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A microscopic approach is presented for calculating general properties of interacting Brownian 
particles under steady shearing. We start from exact expressions for shear-dependent steady-state 
averages, such as correlation and structure functions, in the form of generalized Green-Kubo re- 
lations. To these we apply approximations inspired by the mode coupling theory (MCT) for the 
quiescent system, accessing steady-state properties by integration through the transient dynamics 
after startup of steady shear. Exact equations of motion, with memory effects, for the required 
transient density correlation functions are derived next; these can also be approximated within an 
MCT-like approach. This results in closed equations for the nonequilibrium stationary state of 
sheared dense colloidal dispersions, with the equilibrium structure factor of the unsheared system 
as the only input. In three dimensions, these equations currently require further approximation 
prior to numerical solution. However, some universal aspects can be analyzed exactly, including the 
discontinuous onset of a yield stress at the ideal glass transition predicted by MCT. Using these 
methods we additionally discuss the distorted microstructure of a sheared hard-sphere colloid near 
the glass transition, and consider how this relates to the shear stress. Time-dependent fluctuations 
around the stationary state are then approximated, and compared to data from experiment and sim- 
ulation; the correlators for yielding glassy states obey a 'time-shear-superposition' principle. The 
work presented here fully develops an approach first outlined previously (M. Fuchs and M. E. Gates, 
Phys. Rev. Lett. 89, 248304, (2002)), while incorporating a significant technical change from that 
work in the choice of mode coupling approximation used, whose advantages are discussed. 

PACS numbers: 82.70.Dd, 83.60.Df, 83.50.Ax, 64.70.Pf, 83.10.-y 
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I. INTRODUCTION 

Colloidal dispersions represent one of the simplest classes of materials for which the interplay between viscoelas- 
ticity and externally controlled flow can be investigated. Quiescent dispersions consisting of colloidal, slightly poly- 
disperse (near-)hard spheres exhibit all the hallmarks of a glass transition as the volume fraction is increased. At 
densities above this transition, Brownian motion of the colloids is ineffective in relaxing structural correlations on 
observable timescales: the system remains amorphous, like the fluid phase, but bec omes nonergod ic. The colloidal 



glass transition has been studied in detail by dynami c light scattering measurements (iBartsch et aTl 20021 Beck et all 
I999I, lEckert and Bartsch"2003', 'Hcbrau d et all 1 19971 . iPusev arid van Megen et aLlll987l Ivan Megen and Pusevlll99ll 



van Megcn_a nd Underwood 1993, 1994 ), confocal mi croscopy (IWeeks et alj|2000l) and both linear (iMason and Weita 
W95, Za ckrisson et al.l2006) and nonUnea r rheology ( Besseling et all 20071 Crassous et aLll2006ll2008 . Petekidis et al 
2002, . ,20031 120041 . iPham et aLll2"006l I2008D . The nonlinear rheology of colloidal glasses appears, at least macroscopi- 



cally, to be characterized by the appearance of dynamic yield-stress behavior in which (a) any finite steady shear rate 
restores ergodicity, and (b) a finite limiting stress is attained on slowly reducing the shear rate towards zero. 

Although there is some consensus about the mechanisms (normally described in terms of cage-formation) that 
control the colloidal glass transition, it is not clear a priori whether these also control the nonlinear rheology of 
dense suspensions . For example, many nonlin ear effects have been attributed to ordering or layering of the particles 
(iLaun et al.lll992f ) and/or cluster formation ( Bender and Wagneilll996l iBesseling et al.lT2007l iGanapathv and SooJ 
I2OO6D . Either could be important, especially under conditions where hydrodynamic interactions dominate. This 
dominance seems relatively unlikely at the infinitesimal shear rates seen in colloidal glasses just beyond yield (at least, 
not if the flow remains homogeneous, which we assume here). Several theoretical studies have suggested a connection 
between steady-state nonlinear rheology and the glass transition, but apart from our own work, few of these theories 
explicitly address colloids. For instance, the mean field appro ach to spin glasses was generalized to systems with 
broken detailed balance in order to model 'flow curves' of glasses (jBerthier et al 1 l2000llBerthier and Barr"atll2002D : but 
although the microscopic model is clear, the relation to actual shear flow proceeds only by analogy. The sof t glassy 
rheology model, which describes mechanical deformations and ageing (jFielding et aLll2000l ISoiiich et al. l ll997llSollichl 
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Il998l ). explicitly addresses shear but contains almost no structural information on the material under study. 

Conversely, in a co l loidal context, rheological models have grown up largely without reference to the glass transition. 
For example, iBradvl ( 1993h worked out a scaling description of the rheology of colloids based on the concept that 
the structural relaxation arrests at random close packing (RCP). This contains important insights, but if colloidal 
arrest is actually at the glass transition (which for hard sphere colloids is observed at volume fractions of about 58%, 
well below RCP at 64%), it gives in the end a potentially misleading picture. This is because almost everything one 
could measure diverges at RCP (osmotic pressure, local lubrication resistance, shear modulus...) whereas at a glass 
transition there is only one divergent quantity: the structural relaxation time. This discrepancy therefore cannot be 
fixed by a simple rescaling of the volume fraction to make the RCP point coincide with the observed arrest density. 

In contrast, the mode coupling theory (MCT) of the glass transition gives broadly the right sort of divergence 
in relaxation times (as probed by scattering experiments), without unwanted singularities in pressure or short-time 
dynamics, but does so at a density that is too low by several percent. The results are traditionally shifted to recover the 
correct arrest point, whereafter MCT provides a near-quantitative explanation of numerous dynamic measurements 
at the glass transition in quiescent colloidal dispersions. Two important effects are however neglected; one is aging 



( Purnomo et al.l 



2006^■ the o ther, residual activa ted d ecay pr ocesses at ultra-long times that may cause any glass to 
Gotzelll99ll . rG5tze and Si5grenlll992l[G5Qll999D . Neglect of activated rearrangements defines the 
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so-called 'ideal' glass transition; it is this ideal limit that standard MCT addresses and that wc wish to extend here 
to the case of steady shear. The ultimate flow regime caused by activated hopping would presumably convert the 
yield stress that we find below for the ideal case into a regime of extremely high but finite viscosity (caused by a very 
large finite relaxation time). So long as colloid experiments do not access these ultimate time scales, the yield-stress 
scenario that we develop here should be a good description. 

An appealing aspect of the standard MCT is that its sole input is the static structure factor or equal-time density 
correlator Sq in the equilibrium state (see below). This structure factor is used to create an approximate expression 
for the thermodynamic forces that arise when particles adopt a given Fourier-space density pattern. A separate 
assumption of our approach is that, at a time in the distant past prior to onset of shearing, the system was governed 
by the Boltzmann distribution. This assumption is of course correct for ergodic fluids; however, in the glass it 
represents only one of many possible protocols governing the preparation of the system. Because the system starts in 
equilibrium, slow aging processes (which in practice lead towards that state from a nonequilibrium initial state) are 
precluded from our description. In the context of the present work, which addresses steady-state shear only, this is 
almost certainly not important since, as wc shall see, steady shear restores ergodicity. Thus we expect no dependence 
on the choice of initial state. However, to check this with certainty, calculations would have to be performed with a 
different initial ensemble. MCT-like techni ques t o add ress aging by this route (without shear) have been developed, 
but encounter notable technical difficulties iLatzl (1200X1) . and we do not pursue them here. 

Alongside the present work, se veral MCT-inspired appr o ximat ions have been used to describe the nonlinear 
rheology of colloidal dispersions. Ilndrani and Ramaswamy (Il995l) described self-d i ffusio n at low densities, where 
they suggested a non-selfconsistent, perturbative solution. Mivazaki and ReichmanI ( 20021 ) made this approach self- 
consistent, extended it to collective density fluctuations in dense fluids close to, but always below, the glass transition. 
( Mivazaki et al.l ( 2004 [200(3. ) also presented a field theoretic derivation, evaluated their equations quantitatively, and 
tested the results by computer simulations and experiments.) These approaches investigated the time-dependent 
fluctuations around the stationary state under shear, and in principle required as input the distorted structure factor 
(albeit then approximating this by the undistorted one). In spirit, they followed closely the original MCT without 
shear. However, these authors' wariness of addressing the rheology of the glass phase itself was warranted; their theory 
invokes a fl uctuation-dissipation th e orem t hat cannot be r e hed u pon in the glass. Finally, in the interesting recent 
approach of iKobelev and SchweizeJ ( 2005| ) . ISaltzman et al.l ( 2008t ). entropic barrier hopping prevents glass formation 
in an 'extended MCT' framework, and applied stresses modify the barrier heights. 

In assessing this progress, it is worth recalling some of the prior history of MCT. This has, broadly sp e aking , 
been used to deduce the dynamics from equilibrium structural information in three situations. iKawasa M (119701) 
considered p h ase transitions, wh e re cri tical fluctuations lead to self-similar scaling laws in the structure functions. 
iGotzd ( 199l[ ). ICotze and Siogreiil ( 1992t ) developed MCT for glasses, where the equilibrium structure varies smoothly 
but a bifurcation arises in the equations of motion. Long time tails and back flow phenomena could also be described 
by rather si milar equations, where, however, in contrast to glasses, only hydrodynamic long wavelength fluctuations are 
important (jKawasaki and Gunton.f 973.) . A central tenet for work in these three areas has been that the equilibrium 
structural information used in the mode coupling equations should be under control and well understood. For the 
nonlinear steady-state rheology of a dense fluid close to arrest into an amorphous solid, the equilibrium structure 
structure factor 5*^ is, by deflnition, unchanged from the quiescent case and thus under control. But the same does 
not hold for the steady-state equal-time correlator in the flowing system or 'distorted structure factor'. At a small 
fixed flow rate, the latter quantity, just like the shear stress (to which it is closely related) can be qualitatively different 
just outside and just within the glass phase. Indeed, if the glass transition is accompanied by the abrupt onset of a 
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finite yield stress, both the stress and the correlator should have a discontinuity, in the limit of infinitesimal flows, at 
the glass transition. We shall find that this is indeed a prediction of our approach. (The result is not obvious; for 
instance the SGR model (jSoUich et al.lll997h has a yield stress that rises smoothly from zero on entering the glass, 
giving no such discontinuity.) Thus, for any theory of the nonlinear rheology of the glass transition, a central issue is 
the handling of the stationary structural correlations. We therefore devote part of this paper (Section|V| to discussing 
these in detail. 

Our integration through transient (ITT) approach, first suggested bv iFuchs and CatesI (l2002f ). takes a somewhat 
different route to the challenge of extending MCT to address the rheology of dense fluids close to the glass transition. 
The method proceeds by studying the build up of structural correlations under the combined influence of flow and 
Brownian motion, after switch-on of steady shearing. The initial state, even in the glass, is taken to be the equilibrium 
one, governed by the Boltzmann distribution. (This is also how standard MCT proceeds - for quiescent systems, the 
ideal glass transition is then defined by the loss of ergodicity within an initial Boltzmann state.) ITT differs from all 
the aforementioned approaches to MCT under shear by focussing on the transient density correlators, which follow 
from equations of motion containing the equilibrium structure factor Sq as input. The equilibrium Sq is determined 
by the interaction pair potential among the particles, and is assumed to vary smoothly with thermodynamic control 
parameters. Within quiescent-state MCT, Sq is in fact used as proxy for the pair potential, to calculate thermodynamic 
forces on the particles. ITT follows this avenue directly, rather than making any parallel assumption about the 
nonequilibrium equal-time correlator (the distorted structure factor). This is physically advisable since the relation 
between structure factor and particle interactions only holds in the absence of shear: the structure factor is distorted 
because the system is driven away from Boltzmann equilibrium, not because it adopts the Boltzmann distribution of 
some distorted pair potential. 

The distorted structure factor is then not an input but an output of the ITT approach, just like the stress (to 
which it is closely related). Importantly, we will find - as anticipated above - that the distorted structure factor 
is nonanalytic in density at the glass transition for all finite shear rates, and nonanalytic in shear rate throughout 
the glass. Within ITT, time-dependent fluctuations around the stationary state can be computed, but, in contrast 
with the other approaches, do not play a central role. Specifically their time integrals do not give the transport 
coefficients, such as viscosity, that characterize the steady state. In fact, stationary correlation functions and dynamic 
susceptibilities can be connected via extended forms of the fluctuation-dissipation theorem, but the familiar and useful 
versions t hat apply to linear response around equilibrium states arc violated. (These violations have recently been 
studied bv iKriiger and Fuchs (200S) using approximations beyond those outlined here.) 

Having surveyed the relation to other approaches in the literature, we now summarize the relationship between 
the work of this articl e and our own previous publications on this topic. The first of these was a short paper 
( Fuchs and CatesI I2OO2D outlining in preliminary form: (i) the route via Green-Kubo formulae to exact equations 
that form the basis of the ITT approach; (ii) the use of projection methods on these to obtain MCT-like approxima- 
tions to them in the form of closed equations for transient correlators; (iii) the resulting bifurcation structure; (iv) 
development of semi-schematic and fully schematic MCT models inspired by the closed equations (which otherwise 
remain intractable in three dimensions); and (v) numerical results from these schematic models. The semi-schematic 

(ISHSM) and fully schematic (F^i^) models were subsequently elaborated bv lFuchs and Cate^ (|2003[) w i th additiona l 
results, variants a nd exp e rimental comparisons appearing in several sub s equen t papers: Crassous et al.l ( 20061 . l2008f ) , 
iFuchs and CatesI (|2003bf) . iFuchs and Ballaufj (|2004D . iHainal and Fuchsl (|2008l ) . iHenrich et all (|2005l ). In the present 
work we invoke an ISHSM model in Section |V| (when addressing quantitatively the physics of yielding and of the 
distorted structure factor) which differs marginally from the previously published version, as detailed in Appendix 
C. But otherwise, we do not rehearse any material relatin g to items (iv, v) above. (For comple teness in Section |V| 
we do however briefly restate the results of ( Fuchs and Cates.200^ and ( Fuchs and Catej|2003[ ) on the bifurcation 
structure, item (iii).) 

The major goal of this paper is to give a full exposition of both the ITT formalism prior to its approximation by MCT, 
and the MCT approximations proper: that is, items (i) and (ii) in the above list. However, in developing our MCT- 
bas ed approximations, we make a significant tec hnical change to the ones used originally bv lFuchs and Catei ( 20021 ) 
and lCates et al.l ( 20041 ) . [Fuchs and CatesI ( 20051) . Within the present s cheme, which involves a different definition of 
the transient correlator $q(i) from the one in ( Fuchs and Catej|2002[ ). the initial decay rate Tqit) in the correlator 
memory equation (Eq. (j69p below) is guaranteed positive, not only in the quiescent state, but also under shear. It 
seems desirable to retain this property, since the MCT approach was originally developed under conditions where 
Tq > 0, and both its physical adequacy and its numerical stability are unproven for other cases. Moreover, the 
resulting formulae are generally simpler and more elegant, at least if we neglect certain additional terms which now 
arise in the memory equation, whose form we discuss. Both schemes reduce to standard MCT in the absence of flow, so 
we are free to make this revised choice of approximation. An intriguing consequence of doing so is that it increases the 
mathe mati cal similarity be t ween our mo de-coupling vertex under shear and that proposed bv lMivazaki and ReichmanI 
( 2OO2I ) and lMivazaki et all ( 2004 [2006I ). Given the very different precepts of the two approaches (as detailed above) 
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this does not however imply any deeper equivalence of their theory and ours. 

Our altered choice of MCT closure brings about some mathematical and notational changes which have only a 
minor effect on t h e bas ic structure of the theory, but are pervasive and sometimes subtle. For this reason, although 
iFuchs and Gates! ( 2005[ ) already presented the exact stages of the ITT formalism (item (i)) using the previously 
defined correlator, we re-work much of this material here with the new definition . This has the advantage of making 
the current paper more self-contained, although we still refer to lFuchs and Cate^ (200^ for some important technical 
results that do not depend on the choice of definition made. 

The work presented here on steady states u nderpins two recent short papers in which wc annou nce extensions of ou r 
MCT approach to deal with nonsteady shear ( Brader et al.ll2007) and to general uns teady flows ( Brader et al.ll2008[ ). 
(The first of these uses the original correlator definition bv iFuchs and Cate^ I 2002t ): the second uses the definition 
adopted here.) These rest even more heavily on the ITT approach than does the present work, and are made possible 
because the integration through transients need not assume constant flow rate, nor need the integration continue to 
infinite times. Both simplifications are however retained in this paper which concerns only the steady state, long after 
shear startup. The notational overhead of presenting the ITT method for time-dependent flows is considerable, and 
we have ourselves found the full theory to be much easier to understand once the steady-state version is mastered. 
Therefore we restrict attention to steady shear in this work. 

In common with most MCT-based approaches, we entirely neglect the hydrodynamic interactions that stem from the 
presence of an incompressible solvent surrounding our Brownian colloidal particles. A partial justification for this is the 
hope that, close to a glass transition, the main effect of hydrodynamics is to renormalize the timescale of local diffusive 
transport (as characterized by the bare diffusion constant, Dq). Since at the transition the structural relaxation time 
is a divergent multiple of this local time (as a result of the increasing difficulty, and ultimate failure, in escaping from 
local cages) any smooth density dependence of Dq caused by hydrodynamic interactions is probably unimportant. 
(Note that in the quiescent state, hydrodynamic interactions also cause the Brownian motion of individual particles to 
become correlated, but this does not change the argument.) Under flow, solvent incompressibility also requires locally 
large (but zero-mean) deviations in velocity from that imposed macroscopically. Such deviations could quantitatively 
influence all our results but, for low enough shear rates, need not harm the qualitative picture that emerges. This 
describes the yielding and shear-melting of the amorphous solid or glass. It does not describe any form of shear- 
thickening - a phenomenon which certainly can arise in many dense colloidal suspensions, albeit primarily at high 
bare Peclet number rather than the small ones considered here. Thus our theory should be viewed throughout as a 
low-shear rate approximation; however its validity is not limited to the range of linear response. (Indeed, this range 
shrinks to zero at the glass transition and remains there throughout the glass itself.) While hydrodynamic interactions 
are clearly implicated in some forms of shear thickening, in others its role is less clear. An example of the latter is 
the f ormation under shear of stable arrested granules which then continue to exist on cessation of flow ( Gates et al.l 
I2005D . Elsewhere we have discussed modifications of MCT that can capture these non-hydrodynamic forms of shear 
thickening, but these remain rather ad-hoc and we do not pursue them here (|Holmes et al.ll2003l[2005l) . 

More dangerous is the possibility of macroscopic inhomogeneities in flow rate, for instanc e to form coexisting 
layers of glassy and fluid material, at equal stress but possibly at slightly different densities ( Ballesta et al.l [20081 
iBender and Wagneij|l996llBesseling et al.ll2007llGanapathv and Soodll2006l IVarnik et al.ll2004t) . However, such phe- 
nomena also arise in other fluids such as wormlike micelles (iGates and Fielding 2006f ). and in these cases modelling 
proceeds by first assuming a uniform flow and then analysing the resulting continuum rheology for potential flow insta- 
bilities. This justifies the appr o ach ta ken here which address e s hom ogeneous shearing only. The resulting flow curves 
(explored by fFuchs and Gated (|2003f ) and iHainal and Fuchd (120081)1 generally rema i n mo notonic unless deliberately 
altered (e.g., to account phenomenologically for shear thickening ( Holmes et al.|[2003il2005[ )). This monotonicity rules 
out the most obvious source of shear banding instabili ties, but does not preclud e those involving either coupling to 
concentration gradients, or intrinsically unsteady flow ( Gates and Fielding 20061). Bo th avenues merit further study, 
particularly in view of the recent experimental observations by Ballesta et al.l (|2008l ). which do suggest macroscopic 
flow inhomogeneity under steady shearing in dense colloidal suspensions. 

The rest of this paper is organized as follows. Section [H] details our microscopic starting point and Section Hill the 
exact manipulations that lead to the ITT methodology. Section IIVI addresses the transient density correlators and 
the derivation via MGT of their approximate equations of motion. Section fVl gives a discussion that includes our new 
results for the correlator decay and the distorted structure factor. Section IVTl gives our conclusions; the Appendices 
contain some technical details omitted from the main text. 



II. MICROSCOPIC STARTING POINT 



The system considered consists of N spherical particles (diameter d) dispersed in a volume V of solvent with 
imposed flow profile v(r) = k ■ r, where for simple shear with velocity along the a;-axis and its gradient along the 
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y-axis, the shear rate tensor is k = j xy (that is, Ka)3 = j^ax^fSy). The effect of the shear rate 7 on the particle 
dynamic s is measured by th e Peclet number, Peg = ■^(P /Dq, formed with the (bare) diffusion coefficient Dq of a single 
particle ( Russel et al.lllQSOt ). Dimensionless quantities arc obtained by using d as unit of length, cP/Dq as unit of 



time, and ksT as unit of energy, whereupon Pen = 7. The evolution of the distribution function 5'(r) of the particle 
positions, ri, i = 1, . . . , N (abbreviated into F = {fi})' under internal forces = —diU{T) (with the total interaction 
potential U) and shearin g, but neglecting hydrodynamic interactions, is given by the Smoluchowski equation ( DhontI 
119961 . iRussel et al.lll989D : 

dt^{r,t) = n{T) *(r,t) , 

n ^ n, + 5n = J2^^■{^,~F,~K■ r,) ; (1) 

i 

here fie = J2i ' i^i ~ '^i) abbreviates the Smoluchowski Operator (SO) without shear. In the following, operators 
act on everything to the right, if not marked differently by bracketing. Neglect of hydrodynamic interactions implies 
that we are considering a set of Brownian particles, each of which has mobility /i = Do/ksT; a particle at r feels a 
'flow force' v(r)//i, which for isolated particles exactly replicates the effect of advection, in addition to the interaction 
force Fi from other particles. 

There exist two special time-independent distribution functions, the equilibrium one, and the stationary one, 
which satisfy respectively: 

*e = , f! = . (2) 

The equilibrium one is determined from the total internal interaction energy U via the Boltzmann weight, ^'e(r) cx 
g-£^(r)^ as seen from the useful relation: di ^'e = F.; ^Pe- The stationary distribution function is, however, unknown. 
Equilibrium averages with 'fe will be abbreviated by (...) = J 5'e(r) . . . dV, while determines steady state averages, 
denoted by (. . .)(^) = / ^'^(r) . . . dT. 

The adjoint of the SO can be found from partial integrations (using the incompressibility condition, TracejK;} — 0) 
as: 

fi^ = + + r, • /«^) ■d, = nl + sn^ , (3) 

i 

where boundary contributions will be neglected throughout. We will generally be concerned with the thermodynamic 
limit of 1^ ^ 00 at fixed particle density n = N /V , but where boundary conditions are required, we assume these to 
be periodic. 

The operator 17 acts to the right on a probability density to give the divergence of the resulting probability flux, 
which for f2 = f2e vanishes in the Boltzmann steady state. Accordingly the latter is a right eigenfunction with 
eigenvalue zero. The adjoint operator Sl^ has that same interpretation when acting to the left, but acting to the right 
it represents the flux of a gradient. For = fig there is again a right eigenfunction with eigenvalue zero, but this time 
the eigenfunction is a constant. In practice the transition from a representation involving VL to one involving is 
achieved by partial integration, and is similar to going from a Schrodinger to a Heisenberg representation in quantum 
mechanics. This allows one to work with operators that act on the functions of coordinates whose averages are being 
taken, rather than acting on the probability densities themselves. By this route one obtains equations relating averages 
of different quantities taken within standard (Boltzmann or steady-state) distributions. Such equations then invite 
closure approximations. Closure is instead much harder among equations in which operators alter the distribution 
functions, even if the two formulations contain equivalent information when handled exactly. 

Note further that the manipulations carried out in this Section make no assumption about whether the system is 
glassy " although that does of course inform the choice of approximation made subsequently. Rather, we arc concerned 
here with exploiting the invariance properties of steady states. Time translation invariance restricts the form of static 
and dynamic correlators to those discussed below, whereas translational invariance can also be exploited, as usual, 
by transforming from particle coordinates to Fourier components of the density. The same holds for any spatially 
fluctuating quantity, including the microscopic stress tensor, whose zero-wavevector component is the macroscopic 
stress. 

The shear-dependent operator Jil^ will be shown to capture the afSne distortion of density fluctuations under shear. 
It is given by 

= Y,^^■l^■^^^Y.'^^■l^^ ■ ^^ (4) 
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where Trace{K;} = was again us ed. Without appHed shear the SO fll is a Hermitian operator with respect to 
equiUbrium averaging (lDhontlll996h 



{gnl r) = if* nig) = - Y^i^- ^) , (5) 

i 

and (as seen from speciahsing to f = g) possesses a negative semi-definite spectrum. Here and in the foUowing, 
symbols hke /, g, etc., denote arbitrary functi ons of the full set of particle positions: / = /(F). With shear, however, 
cannot be brought into a Hermitian form ( RiskenlFlQSQf) . This follows from the presence of particle fluxes in the 



steady state. Such fluxes violate time reversal symmetry and there is no detailed balance principle; the latter asserts 
the cancellation of all microscopic fluxes for a quiescent system in steady state. Recently , it has been shown for t he 
present situation, that this entails anomalous fluctuation dissipation ratios at long times (|Kriiger and Fuch^l2008[ ). 

The action of f2 on the equilibrium distribution function fi^e = Sfl'i'e can be written in terms of the stress tensor 
that arises from the interparticle forces: 

,517 *e = - XI *e ~ XI Trace{K; • cr} *e = 7 Cr^y ^e, (6) 

i i 

where the specific form of k for simple shear flows was used in the last equality only. In Eq. is the 

zero-wavevector limit of the potential part of the stress tensor, defined as: 

a^p = -Y,Ftr^^ . (7) 

i 

It is crucial to our approach that the shear stress enters the calculations in two distinct ways: one is the obvious one 
(as a rheological quantity worthy of study) , and the other is as a generator of the transformation between equilibrium 
and nonequilibrium averages. The latter role stems from Eq. ^ and comes to the fore in subsequent developments. 

As detailed in the Introduction, we address only homogeneous (and amorphous) systems so that, by assumption, 
the stationar y distribution fun c tion remains translationally invariant even when it becomes anisotropic as a result 
of shearing. iFuchs and Gates! ( 2005[ ) showed that this assumption is compatible with use of the SO in Eq. ([1]), 



even though the latter appears to break translational invariance. (The proof rests on arguments equivalent to those 
leading to Eq. (|lip below.) At finite shear rate, we shall consider below wavevector-dependent fluctuations around 
the steady state, (5/q = /q — (/q)^^^ and obtain for these quantities not only steady-state averages, but also time- 
dependent correlation Cfg;q{t) and time-independent structure functions Sfg-q (definitions to follow). Translational 
invariance brings appreciable simplifications for all such averages involving time-independent functions of the time- 
varying coordinates F: 

/q(F,t) = e"^*X^/(r)e^'>-"-'. (8) 

i 

Important examples include = 1 which describes density variations (/q(r,t) = pq(t)), and 

from which wc obtain the full wavevector-dependent stress tensor (/q(F, t) — crQ,^(q, {ri{t)})) for particles at positions 
{ri(t)} interacting via a pair potential u(r). While the wavevector depe n dence of density fluctuations is quite familiar, 
the one of stress fluctuations is e.g. discussed bv iBalucani and Zoppil (|l994[ ). As stated previously, the purpose of 
shifting from real space to Fourier space is to exploit translational invariance which means that equal time correlators 
form a diagonal matrix in q space but not in real space. For unequal times, there is a similarly important but subtler 
simplification: correlators can only connect wavevectors that are advected into one another by the intervening flow. 
We next explore in turn these consequences of translational invariance for equal and unequal times. 

Translational invariance in an infinite sheared system dictates that averages involving such quantities are indepen- 
dent of identical shifts of all particle positions: = -I- a for all i, which we denote as F ^ F'. Under such a shift 
the SO becomes 



n\r) =^n\T') -p ■ K- a , withP = X^»- (lo) 

i 
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q(t) 



FIG. 1: Shear advection of a fluctuation with initial wavevector in the x-direction, q(i=0) = q (1, 0, 0)"^, and advected wavevector 
at later time q(f>0) = q (1, —'yt, 0)^ . While Xx is the wavelength in the a;-direction at t — 0, at later time t, the corresponding 
wavelength \y in the (negative) y-direction obeys: A^/Ay — Ax/ Ay — 'yt. At all times, q(t) is perpendicular to the planes of 
constant fluctuation amplitude. Note that the magnitude q{t) = q\/l + isitY increases with time. Brownian motion, neglected 
in this sketch, would smear out the fluctuation. 



Thus, for any fluctuation of a variable which depends on particle separations only, viz. (T) = X/(r'), we have 
PX/(r) = 0. From this it follows that 

/q(r,0 = e-'(''+'i '^')-/q(r',t) , (11) 

whose proof uses the fact that Q. and P • k • a are commuting operators (jFuchs and Catej|2005l ). (This in turn holds 
because the shear rate tensor satisfies k • k = 0, and because the sum of all internal forces vanishes due to Newton's 
third law.) As the integral over phase space must agree whether integration variables F or F' are chosen, it follows 
from Eq. pT|) that steady-state averages can be nonvanishing for zero wavevector only: 

^(/qW><^^=/(7)'5q,0. (12) 

Note that the volume V is taken to be finite at first, with periodic boundary conditions, in order to work with a 
discrete set of wavevectors and with Kronecker J's. Finally, the thermodynamic limit is taken. (For consistency, this 
procedure requires all physical correlations to be short ranged.) Examples of stationary averages are the average 
density n = N/V (which is independent of 7) and the macroscopic shear stress a{j) = {uxy)^^'' /V . 

By similar arguments, wavevector-dependent, anisotropic, steady-state equal-time correlators, built from pairs of 
fluctuations (5/q, (5(7q', are diagonal in q, q' indices. Accordingly we deflne, as one would in a system at rest, 

^/.;q(7) = ^ (-5./* 5g^)^^^ , (13) 

where we adopt the convention that, where no explicit time arguments are given for / or 5, the two times are equal. 
These 'structure functions' arc independent of time in steady state; the familiar equal-time (distorted) structure 
factor, built with density fluctuations, shall be denoted S'q(7) = ((5^?* (Jfiq)^^'. 

The extension of such quantities to unequal times requires explicit account of the time-dependence of the wavevector 
of a fluctuation, as indicated in Eq. (|lip . By translational invariance, the only nonzero averages connect a wavevector q 
at (arbitrary) time t' with its advected counterpart q(t) at later time t'+t (FiglT]). A correlation function characterizing 
this, chosen to closely resemble the corresponding equilibrium quantity, is: 

C/s;q(^,7) = ^('5/qe"'*<Sgq(t))(^). (14) 

Again by convention the suppressed time arguments {t\ say) of /(F) and .g(F) are arbitrary but equal. The exponential 
factor is a time evolution operator, discussed further in Section [llll which defers the evaluation of g, deflned via Eq. 
([8]); until t' + t. (This time evolution operator docs not, of course, act on the q(t) label itself.) 



8 



In Eq. (fTi]) the advected wavevector is defined as 

q(i) = q- q- = q- 7* fey • (15) 



Its magnitude shall be denoted as q{t) = \/q'^ + {qy — '^tq^Y + g^, and its square as q^{t) = [q{t)Y . This describes, 
as in FigUl the advect ion of a density fluctuation through a time interval t. Note that our previous formulation 



( Fuchs and Catesir2002l ) instead utilized the backwards-advected wavevector (q + q • /« t ) . Switching to forward advec- 



tion eases the interpretation and is more suited to the revised mode-coupling approximations that we adopt in this 
paper; see the discussion in Appendix B. (As shown there, this change is not equivalent to merely changing the sign 
of 7.) 

For the special case of density fluctuations, for which /q(r,t) = gq(r, i) = ^?q(^), the abbreviation Cq(t,7) = 

{dg^ e^^ * ^Pq(t))^'^^ shall be used below. This defines the intermediate scattering function of the system under 
shear. Clearly, the complex conjugate relation C^{t,'j) = C_q(f:,7) holds in this case. Because of the inversion 
symmetry of the SO, which, within our approach, is inherited by ^'g and ^'^ by assumption, C^{t,'y) = Cq{t,j) holds 
also. This shows that the stationary intermediate scattering function under shear is real and symmetric in q. 

The steady state averages and correlators defined above all carry explicitly the shear rate 7 as an argument. In 
what follows, wherever this argument is not given explicitly for steady state (as opposed to transient) quantities, 
the shear rate is taken to be zero so that these refer to the equilibrium state, which we continue to assume to be 
homogeneous, and also assume to be isotropic. Thus Sfg-q = ((5/* 6gq)/N is the equilibrium f-g structure function; 
Sq = {Sg^ S0ti)/N is the equilibrium structure factor for the density; and Cq(t) is the equilibrium intermediate 
scattering function. 

The concept of "wavevector advection" , introduced above, is crucial to subsequent developments, and merits careful 
explanation. At one level it is merely a way of subtracting out a trivial effect of flow on density correlations - as 
would be present even for a system with no interaction between particles and no Brownian motion (Fig[T|). Because 
of wavevector advection, if one were to construct a light scattering experiment for such a system under shear, then to 
observe the time correlation of the density (a non-decaying correlation in this example) the detector would have to 
be moved in a specific manner so as to allow for the fact that a particular Fourier-space fluctuation or 'speckle' is not 
stationary but has a deterministic motion in reciprocal space resulting from the shear. (Note that this is strictly a 
thought experiment; light scattering actually measure s an in t ensity autocorrelator from which the density correlator 
follows, in unsheared systems, via the Siegert relation [PusevI ( 199ll) .) Interpreted this way, wavevector advection is a 
book-keeping device for removing this motion. However, below we shall also use the same term to describe the physical 
consequences of advection on the decay of correlators. Specifically, the increase in wavenumber caused by advection 
allows density correlations to be relaxed by much smaller Brownian motions of the colloids than without advection. 
In this second sense, "wavevector advection" represents an actual physical mechanism at work in the shear-induced 
destruction of the colloidal glass. 



III. INTEGRATION THROUGH TRANSIENTS APPROACH 



The expressions collected in Section |TT] for the steady-state properties of dispersions under shear flow require the 
stationary distribution function 'i/g- This satisfies Eq. il'^g = 0; in this equation, time does not enter, and the 
shear rate 7 enters only linearly, via the definition, Eq. ([T]), of the SO. We anticipate however that '^s does not in 
general depend smoothly on control parameters, such as shear rate and density or temperature. In particular the 
ratio {axy)^^y^V = 77(7), computed via '^s, defines the zero-shear viscosity rjQ through the limit rjo = limj^ovH)- 
We expect tjq to diverge on approach to an ideal glass transition, and remain infinite throughout the (ideal) glass 
phase. (If the glass transition is not ideal but rounded by ergodicity-restoring processes, then, nevertheless, a rapid 
variation of 770 takes place.) Thus, in the ideal glass limit addressed in this paper, ^'^ is nonanalytic in 7, and can 
only sensibly be approximated if the glass transition mechanism is taken into account. This mechanism is kinetic, 
not thermodynamic, and consists in the arrest of the structural relaxation at high densities by 'caging' and related 
effects. The ITT formalism for the nonlinear rheology addresses this by calculating '^g via the transient dynamics. 
The time-dependent approach of ^'(i) to the stationary distribution at long times is found from Eq. ([1]), and then 
approximated by considering the slow structural rearrangements of local particle densities, along lines that follow 
closely those developed in the MCT of quiescent glasses. 



A. Generalized Green-Kubo relations 



As detailed in the Introduction, we approach the steady state from a shear startup protocol initialized from a 
Boltzmann equilibrium state at time zero. That is, the system has ^(r,t = 0) = ^'e(r) at times t < 0; at i = 0, a 
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constant shear rate 7 is instantaneously switched on and held constant thereafter: 
The formal solution of Eqs. (|lll6p for t > is: 

^{r,t) = e^'-^^'^eiT) . (17) 

We assume that as t ^ 00, 'I'(r, t) converges to a stationary nonequilibrium distribution ^'s(r), whose solution from 
Eq. p7p can be simplified using the operator identity 

nt _ 1 , / j-L'Mt' 



e'" = 1 + / dt'e'" n . (18) 
Jo 

Combining this with Eq. ([6]) gives for an arbitrary stationary average ( Fuchs and Catej|2002l . [20051 ) 



(7) 



JdT *,(r) /(r) = JdT *e(r) f{v) + ^ j^dt Jdv /(r) e'"^^^' *e(r) a,.y 



= (/>+7 / dtia^ye'''* f) . (19) 







Here, the adjoint SO il^ is introduced via partial integration. The resulting time evolution operator now acts, to the 
right, on the variable / whose nonequilibrium average we want to know, rather than on the probability density. The 
latter becomes time-independent and is therefore given by the initial, Boltzmann, distribution. This formal p rocedure 
is clo sely analagous to the passage from Schrodinger to Heisenberg representation in quantum mechanics ( Messiahl 



Il999f ). In Eq. (|19p . the difference in the stationary and the equilibrium average is determined from integrating up 
the transient correlations between /, the variable of interest, and the flow induced shear stress fluctuation. Taking 
away the test function / leads to the formal expression for the stationary distribution function 

/•OO 

M/, (r) = (r) + 7 dt^,{T)<j,y e^' (r) t (20) 
Jo 

The resulting Eq. (j20p is central to our approach as it connects steady state properties to time integrals formed 
with the shear-dependent dynamics. Knowledge about slow relaxation processes in the system can enter. Moreover, 
averages over the a priori unknown \E's are converted into averages (albeit of more complicated objects) over the 
Boltzmann distribution, rendering them amenable to approximation methods developed for equilibrium dynamics. 

The stationary probability distribution should, like the equilibrium one, be normalized to unity. Therefore we 
require 



0= dt dr ^e(r) a,y ^'>' = / dt dr e''^'>' *e(r) a.,y = / dt {a:,y e'" 1 ) . 
Jo J Jo J Jo 

Here and below, (. . .) (without superscript) represents an equilibrium average; therefore for any constant c 

(a,,ye"'*c)-c(a,.y)=0, (21) 

since the mean shear stress vanishes in equilibrium. The ITT expression of Eq. (|20p is thus confirmed to obey the 
normalization condition. 

Application of Eq. (PD|) , or of the more explicit Eq. , is potentially obstructed by the existence of conservation 
laws, which may cause a zero eigenvalue of the (adjoint) SO, il^ . The time integration in Eqs. (|20ll9p would then 
not converge at long times. This possible obstacle wh en performing memory function integrals, and how to overcome 
it, is familiar from equilibrium Green-Kubo relations (|Forsterlll975 ). One needs to show that the conserved variables. 



which are the eigenfunctions of 17^ with zero eigenvalue, do not cause a non-decaying contribution in the transient 
correlation function in Eq. (|19p . This is achieved by considering the couplings to the densities of the co nserved 



quant ities, which are called 'projections', as introduced by Zwanzig and Mori and others and des cribed bv iForstei 



1975 ). For Brownian particles, the particle number is conserved. Yet, as carefully discussed bv I Fuchs and Gates 



20051 ). density fluctuations do not couple in hnear order to the shear- induced change of the distribution function. 
This follows from the vanishing of the (equilibrium) average {(Jxy e^^^*^ Qq) = 0, and means that no zero eigenvalue 
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arises due to a conservation law ( Fuchs and Catej|2005l ). Because the projection on density fluctuations vanishes, the 



orthogonal or complementary projector Q can now be introduced 

Q = \-P, ^\th P = Y,5b^) ^{5b*^, (22) 

q ' 

which satisfies = Q, P^ ^ P and QP = 0, and in terms of which 

(a,ye"''"*X) = (a^yQe'^'^QX) = {a^yQe'^^'^^^QX) . (23) 

This holds for the observables X = f,g, etc. used in Section[n]to construct steady-state averages, structure functions, 
and correlators. The projection step is exact, and also formally redundant at this stage; but it will prove invaluable 
later on, when approximations are performed. 

From the integration of the distribution function through the transients in Eq. (j20p . we gain explicit expressions 
for the steady-state averages in Eq. (HH) (where, recaU, the definition f{'j) = (/q=o)^^V^ 

was made): 

fii) ^ {U=o)/V + l y^dt (a,,Qe«^^'«*QA/q^o) , (24) 
while corresponding expressions hold for the structure functions from Eq. ()13p . 

SfUi) = (^/q S9^)/N + ^ J^dt {a,y Q * Q A (5/* 6g^) ) , (25) 
and for the fluctuation functions from Eq. p^ . 

Q.;q(i,7) = ('^/qe"'*^3q(t))/iV+^pt'(a,,QeQ"'Q*'QA(<5/*e,"'*55q(t))). (26) 

In Eqs. ((24|) to ([26|) . the symbol AX for the fluctuation in X was introduced, 

AX = X-{X) , thus, e.g., A {Sf* Sg^) = Sf^ Sg^ - NSfg.^ . (27) 

This makes explicit the fact that, owing to Eq. pT|) . all mean values (which are constants, for these purposes) drop 
out of the ITT integrals, leaving only the fluctuating parts to contribute. Note also that all the averages, denoted 
(...), throughout Eqs. (|24|) to ([26|) are evaluated within the (Boltzmann) equilibrium distribution ^'e(r). These 
manipulations may be unfamiliar, but have antecedents in the literature; when studying the nonlinear rheology of 
simple fluids far from any glass transition, generalized Green-Kubo r elations based on transien t correlation functions 
like those in Eg. II24D were found u seful in thermostated simulations ( Morriss and Evang|ll987D and in mode coupling 



approaches ( Kawasaki and GuntonlfTQTSl ). 



B. Transient density correlator 

The problem of calculating steady state averages shall next be converted into one of first finding the transient 
response after startup of steady shear, and then integrating this response in order to use Eqs. ([M)) to To keep 

track of shear-induced structural rearrangements we define the following transient density correlator at wavevector q: 

*q(0 = ^ i^el e"'* (28) 



which differs from the choice made bv lFuchs and Gated (|2002f ) as detailed in Appendix B. Note that ^q{t) depends 



on the shear rate 7; for notational brevity, this argument is suppressed. Its interpretation follows from the joint 
probability W2(Tt,T't') that the system in equilibrium was at state point F' at time t' , and that shear and internal 
dynamics then take the system to state point F at time t > t' . Denoting as P{Tt\T't') the conditional probability of 
evolving via the SO, il, in Eq. ([T]) from state F' to state F, the desired joint probability is given by 

W2iTt, T't') = p{n\r't') *e(r') = e"(^) d{r - r') ^,{r') . 

To derive this, we used the fact that the transition probability P also obeys the Smoluchowski equation ([T]), and that 
its formal solution can be given using the initial condition that F and F' coincide: P{rt\r't) = (5(F — F'). The transient 
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correlator is now defined as the stochastic overlap between an equilibrium density fluctuation Sg^ with wavevector q 
at time t' = and an appropriately shear-advected density fluctuation Sg^j^t^ at later time t, where the time evolution 
between these times is given by the full SO (and so allows for the presence of shear, interactions, and diffusion): 

^''^^^ " A?^ ydr|dr'VK2(R,r'o)<5^,;(r')<5^,q(t)(r) 

Partial integrations, which bring in the adjoint SO fi^, followed by integration over F', then lead directly to Eq. (j28p . 

From its definition and inversion symmetry, it follows that our correlator is real and symmetric in q: <&q(^) = 
$q(i) = ^-q(i)- In the absence of both Brownian motion {Dq = 0) and particle interactions, this correlator does 
not decay at all (3>q(i) = 1); a density fluctuation which is solely advcctcd is tracked perfectly by the wavevector 
advection that was built into the correlator definition, Eq. (|28|) . 

We note in passing that the transition probability P(rt\T't') (which was used not only to obtain Eq. ([28|) but also 
in the corresponding definition of the time-dependent correlators in Eq. (|14p ) is connected to the time-dependent 
distribution function \E'(F,i) via the relation (for t > t' > 0): 

*(r, t)= [ dV p(Tt\r't') ^(r', t') = / dr' e"^^) s{r - F')e"(^') *' ^'e(r') 



gO(r)(i-t') gn(r)t' 



^'.(r) = e"(r)* *e(r) . 



The probability of occupying state point F at time t is thus given by the product of the probability of being at F' 
at earlier time t' , and the transition probability during the time interval t — t' , integrated over F'. Reassuringly, the 
formal manipulations used within ITT conserve the consisten cy between one- an d two-time probabilities required by 
the Chapman-Kolmogorov relations for Markovian processes ( van Kampen|[l98l| ). 



C. Coupling to structural relaxation 

The ITT expressions, Eqs. ((24| to pG]). leave us with the problem of how to approximate time-dependent correlation 

functions of the form {a^y Qe*^^**^* Q AX). Here, as before, AX denotes a general fluctuation. (In the case of the 
temporal correlators Cfg-q{t',j) in Eq. (|26p . AX depends on another time, t'.) The interpretation of these hmctions 
can be learned from their definition 

(cr^y ge'S^^^O* gAX) = jdr jdV W^{TtXQ) a^yiT') AX{T) . (29) 

Here, the joint probability W2 {Tt^T'Q) is formed identically to W2(ri,F'0) in the previous Sect. IIIIB| except that 
the time evolution is given by the reduced SO QflQ, and that the overlap of fluctuations with density at both initial 
and final times is eliminated (again using the projector Q). At time zero, an equilibrium stress fluctuation arises; the 
system then evolves under internal and shear-driven motion until time t, when its correlation with a fluctuation AX is 
determined. Integrating up these contributions for all times since the start of shearing gives the difference between the 
shear-dependent quantities and the equilibrium ones. During the considered time evolution, the projector Q prevents 
linear couplings to the conserved particle density from arising. As stated previously, this projection is optional within 
the current (exact) manipulations since no such couplings arise within the exact dynamics. But in any approximate 
dynamics they might arise (in which case their removal by projection is required, to avoid artefacts). 

The time dependence and magnitudes of the correlations in Eq. ((29|) shall now be approximated by using the 
overlaps of both the stress and A A" fluctuations with appropriately chosen 'relevant slow fluctuations' . For the dense 
colloidal dispersions of interest, the relevant structural rearrangements will be described as usual in terms of density 
fluctuations. Because of the projector Q in Eqs. to ([^5]) . the lowest nonzero order in fluctuation amplitudes, which 
we presume dominant, must then involve pair-products of density fluctuations, PkPp- (These are familiar elements 
in the MOT for the quiescent glass transition.) In accord with the interpretation of Eq. ([29|) . we choose to take a 
static (equal-time) overlap between the fluctuation AX and such density pairs, whose time evolution in relation to 
the earlier stress fluctuation is then approximated using the transient density correlator defined in Eq. (|28p . 
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To present the general case, we give here the mode couphng approximation for the time-dependent transient 
correlator of any two variables, / and g (which may themselves be composite quantities) that do not couple lin- 
early to density, and that evolve in time according to the projected (or 'reduced') SO, Qfi^Q. First, the example 
((5/*(._j^ Q e*^^* Q Sgq) is considered. A projector onto density pairs is introduced as 

■f2 = ^ Sgu Sgp) ^ {Sgl Sg*^ , (30) 

k>p ^ 

where the Gaussian approxima tion, ( Sgt Sg^ (5f3k' Sgp') ~ A^^ Sk Sp (5k, k' Sp.p', was used to simplify the denominator. 
(This is again standard ( Gotzel[l99l[ ).) The ordering k > p and k' > p' should be kept in mind. Fluctuations are 



approximated by their overlap with density pairs, as follows: 



/Xf* ^ Qn^Qtr^, N^.Y^ (^/q(-t) Q ^gk' (5gpO got Qt. . {6Ql5g*p Q Sg^) 

k>pk'>p' P ' ^ 



where the last term vanishes unless q k + p. A crucial step in the mode coupling theory is now to (a) break 
the four-density average into a product of pair averages; and (b) replace the reduced dynamics with the full one 
(iGotze and Siogrenlfl99^ . lKawasaki|[l970l ) : 



{Sgl, 6g*p, e-^ « * Sg^, Sgp) « A^^ S^^.t) *k(-t) (<) %i-t) (t) <5k(-t),k' Sp^.t),p' ■ (31) 

The S'fe(_() are the equilibrium structure factors evaluated with the (magnitude of the) time-dependent wavevector 
k(— t) = k(l + Kt), and capt ure the afhne st retching of equilibrium density fluctuations; see Appendix B for their 
different handling here and bv lFuchs and Gates (,2005i) . GoUecting all terms, and enforcing the wavevector restrictions 
following from translational invariance, we obtain 



QSg^)^ 2^ j^^^ 



— $k(-t)(0 %{-t){t) 



k>p 

k+p=q 



V- ^^f^i-t) Q ^gk' Sgp') (^£'k'(t) Sg*,^f^ QSg^) 

= L c $k'(t)VW- 

k'+p'=q(-^t) 

The last equality follows from a change of dummy-summation indices from k to k' = k(— t) and from p to p' = p(— t); 
the other wavevcctors are shifted from k to k'(t), etc.. A similar shift can be performed in a number of analagous 
summations discussed below, but for brevity this will not be notated explicitly in each case. 

The above mode coupling procedure can be summarized as a rule that applies to all fluctuation products that 
exhibit slow structural relaxations but whose variables cannot couple linearly to the density: 

k>p ^ 

The fluctuating variables are thereby projected onto pair-density fluctuations, and the time-dependence follows from 
that of the transient density correlators ^q{-t)(t)- These describe the relaxation (caused by shear, interactions and 
Brownian motion) of density fluctuations with equilibrium amplitudes. Higher order density averages are factorized 
into products of these correlators, and the reduced dynamics containing the projector Q is replaced by the full 
dynamics. The entire procedure is written in terms of equilibrium averages, which can then be used to compute 
nonequilibrium steady states via the ITT procedure. 

A second rule is needed for fluctuations that can couple linearly to densities. This is derived in complete analogy 
to the case of pair-densities just discussed: 

O+t ^ r \ '^q(-*)(^) / X * fQ1\ 

e « 2^ Ofi-qi-t) ) ( Sg^ (33) 

q 

As a check of Eq. (|33p . we note that it is consistent with the definition of the transient correlator as can be seen from 



NS, ^qit) = {Sg*q e'''' Sg^^,^) « (-^gq ^^^'i-^) ) ^^'ns^^ 



Sg\i Sg^i^t)) 
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= E 

q' 



NS, NS^^,^ V,q(0 = *q(0 



where the Kronecker-i5's lead to the expected result because q'(— t) = q' • (1 + nt) = q{t) ■ (1 + Kt) = q. A similar 
consistency check can be applied to Eq. ((32)) for the pair-product density fluctuations calculated in detail above. 

1. Steady state averages 

The mode coupling approximations introduced above can now be applied to the exact generalized Green-Kubo 
relations from Sect. IIII Al Steady state expectation values from Eq. ([M)) . for variables that do not couple linearly to 
density fluctuations, arc approximated by projection onto pair density modes, giving by the first rule discussed above 

/(7) » (/o>/v + i pt Y. "'t:^^' vi w . (=>") 

with t the time since switch-on of shear. To derive this, the property <I>J^ = <l>_k = 'i'k was used; also the restriction 
k > p when summing over wavevectors was dropped, and a factor ^ introduced, in order to have unrestricted sums 
over k. Within Eq. (|34p we have already substituted the following explicit result for the equilibrium correlation of 
the shear stress with density products: 



{<7xy Q (Jfi-ki-t) <5£'p(-t)) = ^^^^Zf^ ^'k(-t) '5k(-t),-p(-t) - — dtSk(-t) Si^.-p , (35) 



which can be calculated as in equilibrium ( Gotzelfl99l[ ) using '^e and Eq. ([6]), considering time t in the advected 
wavevectors as fixed parameter. The wavevector derivative appears as S'j, = dSk/dk; the second equality, involving the 



time derivative, is useful and will be discussed further in Section |Vl Equation as derived via the mode-coupling 
rule detailed above, contains a 'vertex function' V^, describing the coupling of the desired variable / to density pairs. 
This denotes the following quantity 

= {SglSgu Q Afo)/N = {SglSgu A/o)/iV - ( Sk + n ^) ^) , (36) 



dn J dr ^ ^ 

To obtain the second result, two thermodynamic equalities were used: the first is [Baxter! ijlQG^ l's relation 
(Sgl-Sg^Sgo) = NSq [Sk + n^^), the second is a result for the thermodynamic derivative, A/q)/(|5gqp) 

^q:^)t for '? ~ ^ ' ill which = {,fo)/V is viewed as an (intensive) thermodynamic density (jForsted Il975l 
iGotze and Latj|l989t) . The calculation of the term in Sq is technically somewhat delicate; while this term has 
no effect on the calculation of rheological properties such as stress, it does influence the distorted structure factor 
addressed in the next Section. 

The general result, Eq. (|34p . can now be applied to compute any stationary expectation value, including for example 
the shear stress, {uxyY^^ /V . From Eqs. (j34p to (pH)) one finds the vertex V^'^^ = (k^k„/k) Sl, which give s the explicit 
ITT approximation for the stationary shear stress of a homogeneously sheared dispersion (|Fuchs and Catesii2002i) 

keTj f <Pk klkyky(-t) ^fc^-t) 

/o 7(271)3 kk{-t) Si 



.(,) . . ^ / ^ • (37) 



2. Structure functions 

The structure functions from Eq. (j25p can be approximated along identical lines to the stationary averages and 
become 

r°° k k (—t)S' 

Sf,-Ai) « {Sf^Sg^)/N + ^ jdt ^ ' ' '^'k(-.)(i) , (38) 

with the general vertex built just as in Eq. (j36|l : 



Vl' = {5gl5g,,QA {5f^6g^))/N . (39) 
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Applying this general result to the important case of the distorted structure factor under shear requires use of the 
vertex built with density pairs, which can be evaluated as in Eq. 



V^k^'" = {Sel^l^k Q A {5Q*^Sg^))/N = 2NSl 5^^ - So (^S, + ^ (n5,) . (40) 

The first term results from a Gaussian decoupling approximation of four density fluctuations into products of pairs; its 
factor 2 follows from Wick's theorem (including symmetry), and cancels the factor ^ introduced in Eq. ([M)) in order 
to have unrestricted sums. The second term again follows from Baxter's relation, or alternatively the thermodynamic 
derivative as given after Eq. (|36p . This vertex leads to the ITT approximation: 

S^ii) = S, + ^ y^^^dt <i>^(_t)W} (41) 

where the second term on the right is anisotropic, but the third is isotropic. This result satisfies S'q(7, n ^ 0) ^ 
1 + 0{n) but not S'q^oo(7) — > 1 , as required when only the self-correlations survive in both limits; in Eq. ([4T|) . the 
isotropic contribution does not vanish for q — > oo. Apparently, the mode coupling approximation breaks down at large 
wavevectors where it cannot properly resolve the very local correlations, and misses the fact that all intermolecular 
contri butions to >S'q(7) should vanish there. A previous expression for the isotropic term published bv lHenrich et all 
( 20071 ) inadvertently assumed that this limiting behaviour would be correct; for now, the form above replaces that 



result as the formal prediction from ITT/MCT. As mentioned previously, however, this term stems from the somewhat 
delicate 5*0 piece in Eq. p6p and it is possible that an improved treament will later be found that can restore full 
consistency to this aspect of the theory. Because of the uncertainty over the q range for which this error becomes 
i mport ant, we do not discuss further the isotropic term in the following. However, in numerical solutions, iHenrichI 
(l2007bD : ound it to be subdominant to the first and second term on the right hand side of Eq. (pij) for all wavevectors 



below the second peak in the equilibrium S, 



3. Stationary temporal correlators 



Temporal correlators Cfg-q{t, 7) from Eq. (j26p . describing the stationary, time-dependent fluctuations in the sheared 
system, are next approximated. We refrain from giving the general unwieldy expressions, but note that the derivation 
follows the same method as detailed above, except that a linear coupling of e^** to density is possible in C/g;q(<,7). 
Whenever this coupling does not vanish for specific reasons, the general approximation for the C/g;q(t, 7) follows from: 
(i) using Eq. ((33|) on the i-dependence in Eq. ([26]) . and (m) using Eq. ((32)) on the i'-dependence there. (The former 
is the time delay in the correlator, the latter is the dummy variable in the ITT integral.) As a concrete example, we 
state here the resulting ITT approximation for the stationary density correlator under shear: 



Cq(i,7) « Sq $q(0 + ^ dt' {a,y Q e^^^'Q*' Q5q*^ Sq^.^^^^ 



NSa' 



{5q\, 5gq(t)) 



*q(0 - ^q(7) *q(0 • 



(42) 



It should be emphasized at this point that the stationary correlator Cq(t,7) is conceptually a quite different object 
from the transient one, $q(t) (whose dependence on 7 is, we recall, notationally suppressed). Specifically, the first 
describes correlations between two nonequilibrium states, separated by t in time, with both states arising long after 
startup and in the stationary regime. The second describes correlations between an initial equilibrium state at time 
zero and a transient nonequilibrium one at finite time t after startup. Nonetheless, within our approximation scheme, 
the stationary density correlators differ from the transient ones only by a static renormalization of the amplitude. 
Moreover, this renormalization coincides with the distortion of the structure factor under shear, which is of course the 
initial value Cq(t = 0,7). Therefore, if both are normalized to unity at time zero, the transient and stationary density 
correlators precisely coincide within our chosen MCT approximation. This is arguably a drawback, since discernible 
differences between transient and statio nary correlato r s at i ntermediate times have been observed in both confocal 
microscopy and co mputer simulations by Zausch et al.l ( 20081 ). Very recently, an extended approximation scheme has 
been suggested bv lKriiger and Fuchsl ( 20081 ). which gives a more faithful account of these differences. However that 
extension is not immediately generalizable to other quantities, so we do not pursue it here. 
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IV. TRANSIENT DENSITY FLUCTUATIONS 



The transient density correlators ^q{t) are defined by Eq. ([28]) above, viz., 

where, from now on, in order to simphfy the notation, density fluctuations 6gq wiU be denoted by giq, dropping 
the 6, which is anyway redundant for q 7^ 0. The time evolution operator causes decay by a combination of shear 
and Brownian effects, while the wavevector advection in gq{t) accounts for affine particle motion. Thus the *&q(t) are 

matrix-elements of the time-evolution operator e'^** sandwiched between on the left and Pq(t) on the right side, with 
averaging performed over the equilibrium (Boltzmann) distribution ^ . The usefulness of the ITT approach 
rests on the availability of good approximations to these $q(i), whose exact equations of motion shall be determined 
next. Projection operator manipulations, following standard procedures for systems close to equilibrium, are thereafter 
employed to set up equations of motion suitable for closure by a subsequent mode-coupling approximation. The aim 
is to create an approximate description in which the slow structural dynamics, which is certainly present in quiescent 
situations, can compete with shear-driven relaxation caused by wavevector advection. 



A. Equations of motion 

The time-dependent advection of the wavevector for a density fluctuation can be written using the operator Sil^ 
from Eq. g]): 

gqit)) = e-''''' eq), (43) 

as can be seen easily from 



e 



(The 'kef symbol, ), emphasises that, since the time evolution operator acts on everything to its right, further terms 
will arise if any other variables lie to the right of Qq.) Thus, the operator Sft'^ is identified as causing affine advection 
in density fluctuations. It possesses this interpretation also for higher powers of density fluctuations, e.g. 

^-SnU g^g^'j —"^^g-ityid/dxi gj:(q+k)Ti ^ ^-^t(yid/dxi+yj d/dxj) giq-ri+ik-rj 

= ^g-7ty.(9.+fcx) gi(q-Hk).r, ^ ^-^t(y,q^+V,k^) ^iq-r,+i]^-r, ^ ^^^^^^^^^^y ^^^^ 

Under averaging with the canonical distribution, the advection operator Sil,^ is not self-adjoint. Thus the operator 
g-ao^t gjygg fj^Q advection of a density fluctuation only at the right side of a correlator. Another operator, e''^^**, gives 
the advection of a density fluctuation at the left side: 

(Qqit) = {9*q , with = ^ r, • • (9, + F,) . (45) 

i 

(Note that (5f7t ^ Sil; the latter was defined in Eq. ([T]).) When proving Eq. (PSj) . the presence of the Boltzmann 
weight in the average needs to be recalled, which leads to the differential equation: 

9t (Qqit) = E(*q ■'^ ■'■^ g-'q-(i-*)-r. =-Y^Jdrr,-n^ ■ {d, - F,) ^,*(,)Vi,, = (^^^^^ W , 

i i 

where the last equality follows from partial integration. Its solution is given by Eq. (|^5)) . Analogously to Eq. (|44p . 
the operator also describes the affine motion of higher order density fluctuations at the left side of an average, so 

that(f?;(,) ^'k(t) = (^?q^'ke^*• 
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The transient density correlators of Eq. ((28)) can now be rewritten using Eq. (1431) ^ 

'^^^^^ ^ . (46) 

Introducing the abbreviation U{t) = e^^* e"''^^*, we find by differentiation of Eq. (|^S)) the foUowing result: 

dt U{t) = dt (e"'* e-*'"'*) = e"'* {O^ - 5n^) e'''''' = e'''' 0+ e"*^^** . (47) 

Interestingly, the equilibrium SO appears in this expression. However, because ^ll and 6^l^ do not commute, the 
time dependence of the correlators under shear cannot be written in terms of alone (in complete accord with 
physical expectation). Nonetheless, the appearance here of the Hermitian equilibrium SO suggests we might search 
for another 'well behaved' time evolution operator to simplify the equations of motion for the <I'q(i) prior to subsequent 
approximation. 

As will be argued further below, one candidate for a 'well behaved' time evolution operator is the following: 

ni{t) = e^* nl e-^^'* . (48) 

For density fluctuations and arbitrary products of them, this possesses identical matrix elements to the equilibrium 
SO, ri)j, with the sole difference that all wavevectors are replaced by their time-advected analogues (see e.g. Eq. 
([59)) below). The SO i^l^it) in fact arises naturally in the correlator equations if one cho oses a projection onto time- 
dependent density fluctuations, which is the route we now follow. (This differs from ( Fuchs and Catej|2002t) : see 
Appendix B.) This projection is achieved with the time dependent (Hermitian) operator 



q ^^"^gCt) 

which clearly obeys P(i)^ — P{t). The projector P{t) can be rewritten as 

P{t)=Yl e~'''''^^)^^{s*^^'^'-^e''''''Pe^\ (50) 
where the 'rescaled density projector' 



q 



is introduced. The form of this is very close to P{0) = P, which is the familiar equilibrium projector of Eq. (j22p . 
However, P is normalized differently, as a result of which it is not a true projector (P^ ^ P). 
Using P{t) and its complement Q{t) = 1 — P{t), Eq. (|47)) can be rewritten as 

dt U{t) = e^'* {P{t) + Q{t)) nl e-*'"'* = (e"'* e"^"'*) [PnUt) + nl{t)) = U{t) {PnUt) + nl{t)) , (51) 

where, as promised, the SO n,\{t) now appears. So does a further SO which represents a 'remainder' term: 

nl{t)=e^'''* Q{t)nl e-^^'' . (52) 

This can in turn be separated into two pieces, ilj(t) ~ n,Q{t) + ill^{t), one that does not couple linearly to density 
ffuctuations (Pfig(t) = 0) and another that does (Pr2^(t) 7^ 0). As shown in Appendix A these operators obey 

nl^{t) ^ e^* Q{t) nl e-*'^^'* , nl.{t) ^ e^* S(t) Q{t) nl e-*'"'* , (53) 

ft 

E(t) = 7 / dt' e-*"^*' a.,y e*"^*' . (54) 
Jo 

Note that nglt) reduces to Qnl in the absence of shearing, whereas f^^(i) vanishes altogether in that limit. 
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The decomposition of dtU made in Eq. (|5T|) is thus far purely formal. However, we will find b elow that standard rea- 
sonin g in the framework of projection operator manipulations can be applied to this equation ( Kawasaki and GuntonI 
Il973h . and that it yields, after mode coupling approximations, a numerically stable, self-consistent equation of motion 
for the transient density correlator. (Notably, this holds even though the SO ^\{t) under shear does not live in the 
space perpendicular to linear density fluctuations.) 

We now follow this standard reasoning, whose first step is to introduce a reduced time evolution operator Ur{t,t') 
that satisfies the homogeneous differential equation for t > t' 

dtUr{t,t')^Ur{t,t') nl{t), 

with initial value Ur{t,t) = 1. Its formal solution can be given employing the tim e-ordered exponential functio n e_ 
(see Appendix A) where operators are ordered from left to right as time increases ([Kawasaki and GuntonI [T97I 

Urit,t') = e{''''^^'(^' = exp_ I f'ds (e*"'^ Q{s) ^1 e-^^'^)| . (55) 



Second, by appealing to the uniqueness of the solution of the differential equation ()5ip . and by considering the term 
containing P as an inhomogeneity, it follows that 

u{t) = Ur{t,o) + f dt' u{t') p nl{t') Ur{t,t') . 



Taking a time derivative of this equation leads to the following useful relation for the time evolution operator U{t) 
appearing in Eq. (|46p : 

dtU{t) ^Ur{t,o)nl{t) + u{t)pnl{t)+ f dt' u{t')pni{t')Ur{t,t')nl{t) . (56) 



In a third step, the Zwanzig-Mori-typc equation of motion for the transient density correlator follows from this by 
taking matrix elements between density fluctuations g^. Retaining all non- vanishing terms then gives: 



5t$q(0 + rq(t) $q(t) + f dt' M^(t, t') $q(t') = Aq(t) 

Jo 



(57) 



Note that the term on the right hand side of Eq. ([5S|) does not contribute in situations close to equilibrium, but 
remains nonzero in the present case. This is because the projector Q{t) eliminates coupling to density only at i = 0; 
it does not eliminate the S(t) contribution in Eq. (|53p which subsequently acquires a nonzero value. The resulting 
contribution is denoted in Eq. (j57p as 



Aq(t) = ^ {Q*^Ur{t,Q)^At)Q^)^dt{Q*^Ur{mQ^) , (58) 

which indeed has the property Aq(i — > 0) = (gi* Q Qq) = 0. The presence of this term stems directly from our use here 
of the time-dependent pr ojection P(t), an d repre sents (to us) the only apparent drawback of the current approach 
compared to the one bv lFuchs and Gated (|2002t ). The advantage of the revised approach is however revealed on 
examining the instantaneous friction or 'initial decay rate', rq(i), in the memory equation Eq. (|57p . This is defined 
by the matrix clement of ^a{t): 

(g* nUt) gq) _ {Qq^t) gq(0) _ q'{t) 

Thus, Tq(t) coincides with the equilibrium result ( PusevlflOQll ). where, however, the advected wavevector replaces 
the static one. This expression correctly recovers the 'Taylor dispersion' familiar for non-interacting particles (where 
Mq{t,t') = Aq(i) = holds). Moreover, it guarantees that the initial decay rate in the memory equation Eq. ((57|) 
cannot become negative, since it always corresponds to the initial equilibrium decay rate of some wavevector. (This 
eliminates an important source of numerical instability once mode coupling approximations are applied.) Finally, the 
memory function in Eq. (|59p is given by 

Mq{t, t') = — ^ (g; nl{t') Ur{t, t') nl{t) gq) , (60) 



18 



and corresponds to a generalized diffusion kernel. 

We now further rearrange these exact equations by connecting the generalized diffusion kernel to a generalize d 
friction kernel mq{t,t) which we expect to be amenable to mode coupling approximations ( Fuchs and Catejl2005^ . 
(This step is once again a standard element of the quiescent-state MCT: since both are non-negative, it is easier to 
capture a divergence in the friction than a vanishing of the diffusivity. ) The time-dependent SO is decomposed into 
reducible (fllrit) = nl{t) P{t)) and irreducible (17j(<) = nl{t)Q{t)) contributions, using the non-Hermitian projector 
P{t) defined by 

P{t) = 1 - Q{t) = ^ g^) J {g*^ niit) . (61) 

This is idempotent, P P — P, and exhibits the following couplings to density: P P = P and P P = P, with P from 
Eq. ([22|). (It follows that QP = 0.) Using it, one obtains 

nlit) = Qlit) (Q{t) + Pit)) = Qlit) + ntit) J2 0^) , J,,, , i0*q ^Ut) = ^lit) + ^Ut) , (62) 
where the irreducible part is fll{t) = ^l{t) Q{t) and the reducible one can be written as: 

The decomposition of ^l{t) leads to the following differential equation for the reduced dynamics: 

dtUr{t,t') = Ur{t,t') nlit) - Ur{t,t') (nl{t) + nlit)) , 

which — by the same arguments as given in support of Eq. (|56|) — has the solution 

Urit,t') ^Ui{t,t') + [ dt" Ur{t" ,t')nrr{t")Ui{t,t") (63) 



where the irreducible i Cichocki and Hes"3ll987l . [Puchs and Cated[2005l lKawasaki|[l995l ) fast dynamics Ui{t,t') corre- 
sponds to the solution of the corresponding homogeneous equation. This is given by 



[/,(t, t') = e{*' = exp_ 1^ ds e*'"'^ Q[s) nl e"*'"'* Q(s)| . (64) 
Inserting the expression for Ur{t,t') into the definition of AIq{t,t') in Eq. ([50)) leads to 

A/q(<, t') + rq(i) mq(t, t') r^{t') + Eq (t) Jdt" m^{t , t" ) M^{t" , t') = 0, (65) 
where the friction kernel is defined as: 

"^'i^*' = rmr U'^ "^^^'^ ^r-(t) 0^) ■ (66) 

iV5,(t/)lq(i)lq(t ) 

(Its time-dependence is given by the irreducible dynamics introduced in Eq. (|64p .) 

By an entirely analogous procedure, the term Aq(t) is replaced under the irreducible dynamics by 



Aq(t) = ^ U.it, 0) nlit) g^) , (67) 

which satisfies 



Aq(t) = Aq(t) - Eq(t) / dt' m^{t, t') Aq(t') . (68) 
Jo 

We are nearing the end of our task, which was to write exact equations of motion for the transient density correlators 
in a form suitable for mode coupling approximations. To finish it off, we now view the Zwanzig-Mori type equation of 
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motion ((57)) as a Voltcrra integral equation of second kind for ^q{t), with kernel proportional to Mq{t,t'), and treat 
(Aq(t) — dt *&q(0) /rq(0 an inhomogeneity: 

>^q{t) + dt' ^Mq(t,t') = (Aq(t) - . 

From the theory of Volterra integral equations (jTricomi||l957l ) and the relation of Mq{t,t') to mq{t,t') in Eq. (|65|) . 
one then finds that the solution of this integral equation is obtained with the new memory function mq(<, t'), and is 
given by: 

*qW = (Aq(t) - dt^q{t)) + j\t' mq{t,t') rq(t') ^ (Aq(t') - 9t.$q(i')) ■ 

This result can be brought into familiar form using Eq. (|68p to yield the final, formally exact equation of motion for 
the transient density correlators: 

5t$q(t) + rq(t) |$q(<) + ^ dt' mq{t, t') Of $q(<')| = Aq(t) . (69) 

This is of precisely the form used in quiescent-state MCT, except for the A term on the right. The shear rate enters 
in this term, and in the friction kernel mq(i, i'). Suitable approximations for these quantities are discussed in the 

next Section. If $q(i) is replaced by $q''(t), as defined bv iFuchs and Gated ( 2002f ) and considered in Appendix B, 
and if the irredu c ible d ynamics is obtained from time-independent projections, then the previous ITT equations of 
iFuchs and Gated (|2005l) ^ are recovered. 



B. Mode coupling vertex 

Equation (|69p is exact but can only be evaluated after making approximations for the generalized friction kernel 
(memory kernel) TOq(t,t'), and for the coupling of densities to fluctuating forces, Aq(i). We seek ones that will lead 
to closed equations for these quantities and for <&q(i). 

We start by considering the memory kernel, which is is a correlation function of the fluctuating forces acting on the 
densities, and can be written in a more symmetrical fashion 



The coupling to the stress tensor again arises for arguments analogous to Eq. (|54[) . and here is given by: 



E(i) = 7 / dt' 



whereas the irreducible dynamics is governed, using the decomposition of VL\{t) from Eq. ([53]) . by the evolution 
operator: 

n\{t) = {n^Q{t) + nUt))Q{t) . (71) 

We note that r2g(i) Q{t) P = = Pf2Q(t) Q{t) holds, showing that this time evolution operator finally lives in the 
space perpendicular to density fluctuations. Thus, Eqs. (|67l70l7ip describe the dynamics in the space perpendicular 
to (linear) density fluctuations, but also include couplings to densities via the integrated stress tensors E(<) (via the 
n\.{t) term in Eq. ^) and t{t) in Eq. (|70l) . 

Both E terms vanish initially, but on startup of steady shear increase linearly with accumulated strain, ^t. Ac- 
cordingly, their importance depends on whether or not structural relaxation occurs prior to the accumulation of large 
strains. If the shear-induced decay of the density correlators is completed for small values of 7i, t hen such terms 
can b e neglected in a f irst approximation. This w o uld be consis t ent w ith numerical experiments by iMivazaki et al.l 
(|2004l ). IVarnikI (|2006f) . IVarnik and HenrichI (|2006l ). IZausch et all (|2008f ). which find that the steady state structure 
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factor is only modestly distorted when glasses are sheared just beyond their yield stress. (Equivalently, the yield strain 
is small.) Three furthe r arguments for the neg lect of these terms can be given. First, in the alternative approach 



IS small. j ihree turthe r arguments lor the neg lect ot these terms can be given. I'lrst, m the alternative approach 
originally presented by iFuchs and Cate^ ( 20021 ) . which docs not involve time-dependent projections, these terms do 



not arise under shear. It may yet be possible to establish a projection scheme in which they vanish entirely, while 
maintaining the desirable positivity of the initial decay rates Fq, captured by the present approach. Secondly, the 
form of the A term in Eq. (|69p bear s some resemblance to the 'activated hopping' terms that have been examined in 
several extensions of standard MCT ( Gotze and Siogrenlll987t ). whose goal is to replace the ideal glass transition with 



one in which the relaxation time in the glass phase remains finite (but extremely long) . Arguably it is in keeping with 
a rheological theory of the ideal glass transition to neglect all such terms, not just the familiar static ones. Finally, if 
the yield strain is not small (so that distortions to the structure factor remain substantial as 7 ^ 0), MCT's implicit 
use of a harmonic free energy functional for density fluctuations is called into question; the resulting anharmonicities 
cannot be captured by a theory having the static structure factor Sq as the only input. The S terms could then 
consistently be neglected alongside these anharmonicities. With these considerations in mind, from now on we set 

S(t) ~ ~ t,{t) . (72) 

However, we accept that this approximation goes beyond those that parallel the traditional MCT for quiescent states, 
and might need to be improved upon in future work. 

From this assumption, it follows as shown in Appendix A that Ui{t, t') = Uf^{t, t') which immediately leads to 



A, 



(^^;t/f(^,o)^^^(^)^^q)^o. (73) 



Moreover, on neglect of E, only the components of fi^pq(t) perpendicular to density fluctuations enter on either side 
of the time evolution operator in Eq. (|70p . which means that TOq(i,t') can be written more symmetrically as: 

« ^2(fj%7) ^ (^'^(t') Q(t') e"^*' e^"'* Q{t) nl eq(,)) . (74) 

The quantities (5(t)r2|(3q(t) are the fluctuating forces. Because they do not couple to density fluctuations, in accord 
with our mode coupling precepts, we approximate them by their overlap with pairs of densities. The appropriate 
projection, which is time dependent due to wavevector advection in Eq. (j74p . is then 

P2{t)= V ^k(t)^P(«)^^'^k'(f)gp'(t) _ gk(f)gp(t))(gk(t)gp(t) ^^^^ 

k>p'p>p' <^'kwg;wgk'(t)Vw) ^i^p 

where a Gaussian decoupling of the four-point density fluctuation gives the second form, and the wavevector inequal- 
ities prevent overcounting. 

Using this in Eq. (|74p for the fluctuating forces gives 

« q2^t)q^t') ^ ^^^(*') ^^^^'^'^ '''''' (^^^ 

Finally, we factorise the four-point correlation function with reduced dynamics into the product of correlators with 
full dynamics, just as was discussed in connection with Eq. (|3ip : 

(gk(t') QUn e-^*' C/^(t,t') e*"'* gk'(t) Qp'it)) « S^ir) ^K*') *k(t')(^ - *P(t')(^ - ^K^' ^p p' ■ (77) 

The Kroneckcr-(S's arise from translational invariancc and homogeneity. The vertex functions Vqkp measure the 
overlap of the fluctuating forces with the density pair fluctuations in Eq. (|76[) and can be evaluated from equilibrium 
information. Although th ey are r equire d as functions of the advected wavevectors, their evaluation proceeds exactly 
as in equilibrium MCT in ICotzd (|l99lf ) 



{g* nl Q pk Qp) 

^qkp = KTTTT' ^ ^ ■ i^^(^k +pnCp) 6qM+p , (78) 

where Cq is the equilibrium direct correlation function connected to the structure factor via the Ornstcin-Zcrnike 
equation Sq = 1/(1 — ncq). (As is standard practice the convolu tion approx imation has been used in Eq. (j78p to 
neglect a small contribution from higher order direct correlations (jCotz 



21 



The final expression for the memory function in our mode coupling approximation can then be written: 

- ^'g2(^)^g(^y'^ K,.p(t)F,.p(0 %it')it~t') . (79) 

Here p = q — k, and we have abbreviated Vqkp(t) = Vq{t)k(t)p{t) ■ A change of the integration-variable from k 
to k' — k(t') leads to an alternative expression which may exhibit more clearly the various origins of the time- 
dependences. We write: 

mq(i,0 =mq(t,)(t-t') , (80) 

where the reduced memory function rhqi{T) is evaluated at the time-dependent wavevector q' — q(t')i ^^'^ depends 
on time only via the difference t = t — t' . This quantity is given by: 

- ^E%if^ ^.-'p'(-)^.'-p' 'J'p'M ■ (81) 

Here, p' = q' — k' holds, with q(t) = q'(r) and analogous expre ssions for the other waveyector s. 

Notably, the mode coupling vertex derived from field theory ( Mivazaki and Reichmanll2002l . [Mivazaki et al.ll20d3 ) 



for the time-dependent fluctuations around the stationary state apparently coincides with Eq. ([8T|) . albeit with 
the important difference that the distorted structure factor enters there in place of our static one. (The different 
philosophies behind either approach were discussed in the Introduction.) If the 'external' wavevector q lies in the 
plane perpendicular to the flow direction, q • x = 0, and thus is not advected, q(t') = q' = q, then the memory 
fimction simplifies further, as time enters only via the time difference t = t — t' , Tnq-q^=Q(t,t') = "i-ij^ ,5^ (t) . 

The generalized friction kernel of Eq. ([79|l . which vanishes in the absence of particle interactions, scales like 
for small wavevectors. Whereas one power in q arises from density conservation, the second arises because the total 
force among all particles vanishes via Newton's laws. (The present result for rriqit, t') of course recovers the standard 
MCT expression for 7 = 0, to which similar remarks apply.) Viewed as a function of t, the memory function under 
shear has a positive maximum value at t = t', which describe instantaneous friction. As t' is increased, this value 
decreases under the cumulative effects of shearing between startup and time t' . Shear decreases the correlations by a 
dephasing of the two vertex factors which enter Eq. (j79p . The two factors coincide only at t = t' , creating a squared 
vertex familiar from standard MCT. (At other times, the product is not necessarily positive so that stability of Eqs. 
(|69l79p is not automatic.) This dephasing results from a shift of the (internal) advected wavevectors to higher values, 
thus suppressing the effective interaction potentials c^^oo — *■ and decreasing the friction. Additional decorrelation 
during the time interval t — t' enters via the density correlators; this represents a 'Brownian decay factor' (in the 
sense that these correlators do no t decay without Brownian rn otion; see Section ITV|) . The overall effect of shearing, 
as previously discussed elsewhere ( Fuchs and Catesll2002l [20031) is to cut off memory and thereby fluidize the system. 



V. SUMMARY AND DISCUSSION 



Our combined ITT/MCT approach to the rheology of steadily sheared suspensions consists of the approximated gen- 
eralized Green-Kubo relations summarized in Section UlI Al which introduce the transient density correlator (defined 
in Eq. (j28[) ) describing structural relaxation, and its equation of motion, which results from Eq. (|69p : 

dt%{t) + rq(t) |$q(t) + ^ dt' mqi^t'){t - t') df $q(t')| = • (82) 

The initial decay rate is given in Eq. (|59[) and the mode cou pling approximation for the memory function in Eq. (|8ip . 
Since the preliminary presentation of the ITT approach bv lFuchs and CatesI ( 2002f ). a number of r esults found using 



various formulations and simp lifications, based on this general framework, have been worked out ( Fuchs and Catea 
I2OO3I iHainal and Fuch"il2008f ). All these are underpinned by the more complete presentation offered here, which 
however differs from the original version, as mentioned previously and summarized in Appendix B. These differences 
do not affect the main conclusions from the ITT approach, and it lies beyond our scope to present details of all of such 
predictions here. Below, we summarize some of the most important ones, point to the literature for more detailed 
discussions, and compare some additional predictions with recent experiments and simulations. Because the transient 
density correlators $q(t) determine the steady state properties, we start with them. Let us recall, from Eq. ([42]), that 
within the current mode coupling scheme transient and steady state density correlators essentially coincide, being 
connected via Cq(i,7) = >5'q(7) *&q(i). Although as mentioned in Section [ill CI 3 this is only an approximation, using 
it allows comparison with a wider range of experimental and simulation data. 
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A. Universality aspects 



The ITT/MCT equations require as input the equihbrium structure factor Sq, the average density n, and the shear 
rate 7. The dynamics on short time scales, which is not treated adequately in ITT, is specified only by the bare 
(short time) diffusion coefhcient Dq which enters Tq{t) as a factor (we set Dq — 1). The ITT equations contain a 
bifurcation in the long time limits /q of the transient density correlators. <f>q(t 00) = /q. This nonequilibrium 
transition generalizes the ideal glass transition found in standard MCT, and its predicted position as a function of 
density, temperature etc (entering via Sq) coincide with the standard one: there is no shift caused by 7. Nonetheless, 
all nonergodic states become ergodic under infinitesimal shearing, that is, the glass 'form factors' vanish, /q = for 
7^0. Thus, enforcing even an infinitesimal steady shear flow melts the glass, creating a state in which all density 
fluctuations - including those not directly advected by shear - relax in finite time. This follows simply because shear 

advection cuts off the structural memory and forces fhq{t — > 00) = 0. 

A true noncrgodicity transition, like the ideal glass transition of standard MCT ( G5tzelfl99ll iGotze and SiogrenI 



thus docs not exist in ITT/MCT (unless one sets 7 = from the outset, in which one recovers the standard 
theory). However, at the locus of the MCT glass transition, ITT/MCT predicts a transition from a viscoelastic fluid 
to a yielding glass. In fluid states, steady state averages, correlators, etc. show a linear response regime in shear rate 
amounting to a regular Taylor series in 7 (at least in low orders); here, the theory reduces to MCT for 7 — > 0. In 
glass states, however, averages obtained from the generalized Green-Kubo relations do not possess a linear response 
regime. This follows because the final relaxation of the transient correlators now arises only because of the drive by 
shearing. The integrals in the Green-Kubo relations th ereby attain non-analy t ic dependences on 7. (See Sec tion IV CI 
for the example of the distorted structure factor and ( Fuchs and Catej|2002l l2003l iHainal and Fuchs^l2008^ for the 



shear stress.) A posteriori, this result is not surprising, because, as Maxwell explained (and MCT recovers) glasses 
are solids and should thus exhibit a linear response regime to strain (or deformation gradient) but not to shear rate 
(or velocity gradient). 

Close to the glass transition, the ITT/MCT correlators exhibit a two-step relaxation, comprising a decay onto /q, 
followed by the relaxation to zero. The dynamics close to the p lateau /q factoriscs into time-independent amplitudes 
and a time-dependent scaling function (jPuchs and Catej|2003[ ): 

*q(0 = + hq Git) , (83) 

where the critical glass form factor and the so-called critical amplitude hq are isotropic, not affected by shear, 
follow from the equilibrium structure factor at the tra nsition point, and agree with the corresponding quantities from 
quiescent MCT (lGotzelll99lllG6tze and Siogrenlll992f) . 



The scaling function Q{t) carries all the important dependence on time, shear rate and other control parameters. In 
quiescent fluid states, the flnal relaxation, called the a-proccss of glassy relaxation, is isotropic, and possesses a finite 
internal relaxation time, to be denoted as r. This depends on a 'separation parameter' e which is any combination of 
thermodynamic control parameters that vanishes linearly at the arrest transition; for colloids we choose e = (</>— (j)c)/4>c 
where denotes volume fraction and 0c the critical value of this. The scaling function is initiated by a power law, 
called the von Schweidler law Q{t — > c»,£ < 0,7 = 0) oc —{t/r)^, with exponent b smaller than unity. ITT identifies 
the 'dressed Peclet' or Weissenberg number Pe = 7T as the expansion parameter that controls when the effect of 
shearing the fluid starts to matter, at Pe ^ 1 . This of course depends sensitively on e because the a-relaxation time 
does so. 

In states which are nonergodic glasses at rest, the relaxation time t formally is inflnitc. As already made clear, 
such states are melted by any nonzero shear rate 7, as can be deduced from the nonlinear stability equation for Q{t) : 



c(^) i^tf + A g^t) ^ ^ fdt' Git - 1') g{t') , (84) 







which, along with its initial condition 



g{t ^ 0) ^ (Vio)-" , (85) 
was derived via ITT bv iFuchs and Gated (I2002D . The so-called critical law in Eq. ([55)1 is the asymptotic solution 



for long times right at the bifurcation point £ = = 7, which can be matched onto the short-time dynamics by 
chosing the time scale to- The parameters A and c'-'^-' in Eq. (|84p are determined by the static structure factor at the 
transition point, and, in the case of hard-s phere colloids, take valu es around A w 0.73 and c^^^ ~ 0.65 if one estimates 
Sq from the Percus-Yevick approximation ( Fuchs and Catej|2003l ). The transition point then lies at packing fraction 



3 



UcRh ~ 0.52, and the scaled separation parameter entering Eq. (|84p obeys e = Ce with C ~ 1.3. The 
exponent a in Eg . (1551) depends on the 'exponeiit par ameter' A via A — T{1 — a)^/r(l — 2a) where P denotes the 
Gamma function (|Gotzelll99ll . ICotze and Siogrenlll992l ). 
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The presence of a finite shear rate enforces ergodic decay of all density fluctuations, not just those directly advected 
by shearing. Indeed, the long-time solution of Eq. (|84p for all separation parameters e xhibits a finite-time rela xation 
(or 'yield process') whose initial decay from the plateau is governed by the asymptote ( Fuchs and Catesll2002[ ): 



g{t^^)^-J^—j\^t\ = -t/T^. (86) 

The effect of hydrodynamic interactions (HI) on the local dynamics, were this to be incorporated into our approach, 
could be expected to shift the matching time to in Eq. (|85p . (Note that 7^0 ^ 1 is required for the asymptotic 
solution provided by the factorisation law in Eq. (|83|) to become valid.) However our results for the the shear-induced 
relaxation time ~ I/7 in Eq. ([86|) cannot be changed without wholesale alterations to the structure of the theory, 
in which shear only enters via the affine advection of fluctuation patterns. 

Schematic models, in which all or part of the wave- vector dependence in the theory is suppressed ( Fuchs and Catei 



l2003l . lHainal and Fuchsl[2b08[ ) can capture all aspects of G {t) and therefore provide an important toolbox for exploring 
generic aspects of the transient dynamics in ITT. (For the same reason, such models played an important role in the de- 
velopment of standard MCT for quiescent glasses (|Gotzelll984l . Il99ll ).) In particular the F^j' niodel ( Fuchs and Catej 



in which one considers only a single representative wavevector, has been used extensively to analyse the steady 
state flow curves (7(7) obtained in colloidal dispersions of hard and soft spheres and in computer sim ulations of binary 
supercooled mixtures (jCrassous et al.ll2006l l2008l iFuchs and Catesll20Q3l iHainal and Fuchsll2008h . With rather few 
adjustable parameters, very good global agreement between measured steady state data and fltted curves from the 
F^^'' model could be obtained. The semi-schematic isotropically sheared hard sphere model (ISHSM) was developed 
to wholly incorporate the full isotropic dynamics of quiescent MCT, while simplifying the rather complicated shear 
advection effects presented above for the ITT/MCT theory. This was done by selective use of an 'isotropic average ' 
in wh ich all wavevectors are effectively treated as though pointing in the vorticity direction ( Fuchs and Catejl2002l 



I2OO3D . An ISHSM model, whose results are compared to spatially resolved data from experiments and simulations in 
Section IV Bl is summarized for completeness in Appendix C. 



B. Yielding process 

The shear induced decay of density fluctuations causes the transient correlator $q(t) to fall below the plateau 
<f>q ~ fg that would persist indefinitely for an unsheared glass. Since long time Brownian motion is frozen out, we might 
expect this decay to be purely strain-induced and thus rate independent, in which case a scaling law, $q(t, 7) — > <I'q(i), 
should describe the yielding process. Here, the re scaled time is i = t/r -y oc <|7|. The existence of such a scaling law 
has been called time-shear-superposition principle ( Besseling et al.ll2Q07t) . At the transition, e = 0, and neglecting for 



small accumulated strains 7t ^ 1 the external advection in the memory function. TOq(t')(i — t') w mq{t — t'), the 
equation for the master function $q was found as (|Fuchs and Catesll2003l ) 



$q(i) = mq(t) - -~ / dt'mq(t-t')*q(0 ■ (87) 

dt Jo 

The memory function ■mq{i—t') is the one defined from Eq. ([8T|) . evaluated at e = and with the asymptotic scahng 
form $q replacing the actual correlator. The result Eq. ([57)) highlights the accelerated decay of correlations induced 
by shear. The structural decay depends on the accumulated strain t oc \tj\. 

The yielding process is initiated by Eqs. (|83l86p . <&q(t ^ 0, e = 0) = — hqi, which describes an initially isotropic 
yielding process under shear. The linear initial decay of ^q{i) also suggests little stretching of the final decay. This 
suggests an isotropic exponential as approximation to the yielding master function $q(i') ~ fg exp {—t/rlj, q)}, with 
'''Hi l) — ifq/f^q) ''7- However, the true degree of anisotropy of the process and its non-exponentiality at later rescaled 
times is still unknown. 

Note that all quantities characterising the short-time motion, like the matching time to, have dropped out of Eq. 
([57)1 . Its solutions $q(?), and the stationary averages obtained from them, like the yield stress a~^, thus in ITT depend 
on the equilibrium structure factor only. Hydrodynamic corrections to the bare diffusivity, for example, are predicted 
not to affect the value of the yield stress, which at the transition can be obtained explicitly from 



+ _ kBT /2A-1 f d'k klkyky{-i) S'kS'^^^i^ 
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FIG. 2; Steady state incoherent int ermediate scatte r ing fu nctions ^q{t) as functions of accumulated strain for various shear 
rates 7; the data were obtained by iBessehng et alj (|2007l ) in a colloidal hard sphere dispersion at packing fraction (j> = 0.62 
(at e ~ 0.07) using confocal microscopy; the wavevector points in the vorticity (z) direction and has q — 3.8/ R (at the peak 
o{ Sq). The effective Peclet numbers Peeff = ' iR^j/ Ds are estimated with the short time self diffusion coefficient Ds ~ Do/ 10 
at this concentration from Ivan Megen et al.l (|1998D . ISHSM calculations with separation parameter e — 0.066 at qR = 3.9 
{PY-Sq peaking at qR — 3.5), and for strain parameter jc = 0.033, are compared to the data for the Feoff values labeled. The 
yielding master function at Feoff ~ lies among the experimental data curves which span 0.055 < Pceff < 0.45, but discussion 
of the apparent systematic trend of the experimental data would require ISHSM to better approximate the shape of the final 
relaxation process. 



While data for the transient density correlators remain scarce (|Zausch et al.ll2008[) . detailed measurements of the 
stationary dynamics under shear of a colloidal hard sphere glass have recently been performed bv lBesseling et ah l (l2007D 
using confocal microscopy. Single particle motion was investigated. We can therefore compare theory to experiment 
for the density relaxation in a shear-melted glass at roughly the wavevector inverse to the average particle separation, 
if wc make three further assumptions: (i) that transient and stationary correlation functions agree up to an amplitude 
factor (Eq. (|1H)); (m) that it is valid to compare the measured incoherent density correlators to the coherent ones 
calculated in ITT; and (in) that the ISHSM (Appendix C) is adequate for such a comparison. Figure[2]then shows self- 
intermediate scattering functions measured for wavevectors along the vorticity direction where neither afRne particle 
motion nor wavevector advection appears. The stationary correlators deep in the glass, for shear rates spanning almost 
two decades, are shown as function of accumulated strain jt, to test whether a simple scaling Tj ~ I/7 as predicted by 
Eq. (1871) holds. Small but systematic deviations are apparent which have been interpreted as a power law r-y ~ ^"O-S 
(iBcsseling et al.|[20 07. Saltzm an et al.. .2008.). (ISHSM computations were performed for a nearby wavevector where 
Sq is around unity so that coherent and incoherent correlators may be assumed to be similar as argued by iPusevI 
( 19781 ).) A separation parameter e very close to that of the experiments was taken for the fit, allowing for a small 
deviation so as to match better in amplitude the final relaxation step. The yielding master function from ISHSM can 
be brought into register with the data measured at small effective Peclet numbers Feoff, by using a phenomenological 
'strain rescaling parameter' 7c = 0.033 (see Appendix C for the definition of this). For this parameter one expects 
values of order unity to be able to compensate for the oversimplified treatment of angle averaging in ISHSM; the 
smallness of the fitted value is not yet understood. The effective Feclct number Pe^s = ^R^i/Dg introduced for 
these fits measures the importance of shear relative to the Brownian diffusion time obtained fro m the short time 
self di ffusion coefficient at the relevant volume fraction. (A value Ds/Dq = 0.1 was taken from Ivan Megen et al.l 
(Il998l ).) The shape of the final relaxation process differs subtly between experiment and theory, and therefore we can 
make no definitive comment at this stage on the observed small systematic drift of the rescaled experimental curves. 
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FIG. 3: Steady state incoherent intermediate scattering functions measured in the vorticity direction as functions of 

accumulated strain ■yt for various shear rates 7; data from mol ecular dynam ics simulations of a supercooled binary Lenard- 
Jones mixture below the glass transition at {Tc — T)/Tc ~ 0.3 (|Varnikll2006l ). These collapse onto a yield scaling function at 
long times. The wavevector is g = 3.55/7? (at the peak of Sq). The quiescent curve, shifted to agree with the one at the highest 
7, shows ageing dynamics at longer times outside the plotted window. The apparent yielding master function from simulation 
is compared to the ones calculated in ISHSM for glassy states at or close to the transition (separation parameters e as labeled) 
and at nearby wave vectors (as labeled). ISHSM curves were chosen to match the plateau value fq, while strain parameters 
7c = 0.083 at £ = (solid line) and 7c = 0.116 at e = 10~^ (dashed line) were used. 



(This drift is responsible for the fitted dependence ~ 7 reported by iBesseling et al. I (|2007D .) At the lar ecr 



efi'cctive Peclet numbers, Pcoff > 0.5, for which the short-time and final (shear-induced) relaxation processes move 
closer together, the model gives quite a good account of the 7-dependence. Overall, theory and experiment agree in 
finding a two step relaxation process, where shear has a strong efi'cct on the final structural relaxation, while the short 
time diffusion is not much affected. 

The shape of the fin a l relax ation step in a shear-melted glass can be studied even more closely in recent computer 
simulations by IVarni^ ( 2006( ). where a larger separation of short and long time dynamics could be achieved. In 
these molecular dynamics simulations of an undercoo led binary Lenard - Jones mixture, schematic ITT models give 
a good account of the steady state flow curves, (7(7) (jVarnik and Henri ch 2006|). Figure [3] shows the corresponding 
stationary self intermediate scattering functions for a wavevector near the peak in Sq, oriented along the vorticity 
direction, at shear rates spanning more than four decades. Collapse onto a master function when plotted as function 
of accumulated strain is nicely observed as predicted by Eq. ([57]) . At larger shear rates, the correlators peel away from 
the master function; this resembles the behaviour observed in the confocal experiments in Fig.[2l Using the additional 
assumptions (z) to [iii) given above, the shape of the master function can be fitted using our ISHSM calculations. As 
before, to bring these into register, a strain parameter 7c was introduced, whose smallness again remains unaccounted 
for. After this rescaling, modest but visible differences in the shapes remain: the theoretical master function decays 
more steeply than the one from simulations. 



C. Distorted structure under shear 



From the transient correlators discussed (and, partially, validated) in the previous Section, all stationary averages 
then follow within ITT-MCT approach. While the re sulting flow curves and yield stresses have been extensively 
predicted via both ISHSM and schematic ITT models (|Fuchs and Gated |2003[) . the distorted microstructure has so 
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FIG. 4: (Isotropic) distortion of the structure factor SSq{'y) under shear within ISHSM approximation in glass states. The 
limiting values for 7 — » (Pceff — 0) are shown at the transition (e — 0, black line), and at e = 0.066 (red line) corresponding to 
the confocal data in Fig. [2] at <f) = 0.62; for the latter e, the increase of SSq{'j) with shear rate is shown (blue dot-dashed line); 
strain parameter 7c = 0.033 was chosen. The upper thin dashed lines indicate the equilibrium structure factors Sq from PY 
approximation for these densities (with corresponding colour code). The inset shows the complete distorted structure factors 
Sq{'j) a.t (fi — 0.62 (red solid line for 7 ^ 0; blue dashed-dotted line at Pecff = 0.6) compared to the equilibrium ones (black 
dashed line); Spline interpolation through the data on the grid of the main figure was used. At e = the difference for 7 7^ 
can not be resolved. 
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far only been calculated within ITT to lowest order in shear rate ( Henrich et al.l[2b07t ) 



As discussed in Section U the microstructure possesses a linear response to 7 only in fluid states, while in yielding 
glass states, the limit S'q(7 — > 0) does not agree with the equilibrium structure factor denoted Sq. Rather ITT-MCT 
provides the general results 

with 

5S-f^°{i) = \^dt (j^Sq^^t)) $^(-t)W (88) 







where the substitution ^Sqi^^t) = 7 '''^qf^'lt)^ '^qi-t) made in order to tidy up the expressions given in Eq. (j4ip . 
As discussed in Sect. IIII C 21 we refrain from discussing 55'^''° (7), which is quantitatively smaller for not too large 
wavcvectors. The rewriting in Eq. (|88p highlights how the distorted structure arises from the affinc stretching 
of (equilibrium-amplitude) density fluctuations, in competition with structural rearrangements as encoded in the 
transient correlators. Equation ([88]) also shows that whenever the decay is shear- induced (so that the final decay of 
<i>q(t) depends on time via t/rj ~ t|7|) the limit SSqi'j — > 0) follows. This hmit can be studied by replacing $q(t) 
in Eq. ((88)) with $q(f) from Eq. fST]) . and by performing the time integrals over the rescaled time i oc t\'^\. For the 
anisotropic contribution for example one then finds the 'yield value' of the structural distortion 

SSf^°-+{e > 0) ^ 6Sf^°{j ^ 0, £ > 0) Jlli *q(-t)(*~) ' (89) 

where the right hand side is explicitly shear-rate independent. The discontinuity represented by dS^{e > 0) 7^ 
reflects the same physics as the dynamic yield stress; we therefore call it the 'dynamic yield contribution' to the 
structure factor. Figure 2] shows distorted microstructures under shear for glass states at infinitesimal and finite shear 
rates. These results were obtained within the ISHSM of Appendix C, at parameter settings corresponding to those 
fit to the confocal microscopy data in Fig. [21 so as to assess the magnitude of the effects in a realistic case. Thanks to 
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the isotropisation within ISHSM, the distorted microstructure becomes isotropic in the yielding glass. Right at the 
transition, e = 0, only a tiny discontinuity SS^ appears. This suggests that the dynamic yield contribution to the 
structure factor might not be discernible until one moves deeper into the glass. 

The connection between the dynamic yield contribution in the stationary structure factor >5'q(7) and the actual 
yield stress becomes clearer from recognizing a mathematical relation between Eq. p7|) and Eq. (|4T|) . whereby the 
shear stress can be obtained from an integral over the distorted microstructure: 



cr(7) 



'^^^ Cfc ky 

(2^ ^'-^^^ ' 



(90) 



where c'j, = dc^/dk enters. The direct correlation function c(r) plays here the role of the negative of an effective 
potential in units of thermal energy. Because of the angular dependence in Eq. (|90p , only the anisotropic contribution 
jj^aniso Qj^i^Qj-g [j^ tjjg calculation of shear stress, giving a direct connection between the yield stress contribution 

to the distorted structure factor (5S'q'"''°'+(e > 0) and the dynamic yield stress itself. (The latter is defined as 
-(£ > 0) = cr(7 ^ 0,e> 0).) 



The relation in Eq. (^0]) was previously found by iFredrickson and LarsonI (1983) in a Gaussian field theory for 
block copolymers. We interpret this similarity as an indication that the ITT approach based on MCT implicitly 
assumes a Gaussian distribution for the time-dependent density fluctuations. This is connected with the breaking of 
higher order density moments into products of time-depende nt density correlation f unctio ns, which lies at the heart 
of our approach. In this context it is interesting to note that ICrooks and Chandleil ( 1997t ) found equilibrium density 
fluctuations in a dense fluid simulation to follow a Gaussian distribution rather closely even for the highly nonlinear 
interaction given by the hard core excluded volume constraint. It would be very interesting to see this tested under 
flow also. 

Let us turn finally to the ITT-MCT predictions for the steady shear stress, which defines the flow curve a{'y); 
the shape of this curve depends solely on Sq (with scale factors in time and stress deriving from Dq, ksT, and 
the particle diameter d). Numerous s uch flow curves have by now been obtaine d, fo r a variety of sy stems with 
differing Sq, either within the ISHSM (|Brader et al.ll2008l . iFuchs and C"ates 2002 . 200l iMivazaki et "alT[2b04.. .20061 



Zausch et al. 2008) or via fu l ly schematic models like the F^2^-model ( Crassous et al.l [20061 l2008l iFuchs and Cate^ 
20031 iHainal and Fuch^ '2008'. Varnik and Henrich '20 06|) . Note that exactly the same ISHSM as used here (Appendix 



C) forms the basis of Refs. (.Brader et al.. .2008. Zausch et al. 2008); this differs in de tail from the earlier version, as 
used in Refs. (|Brader et al.ll2007l . iFuchs and Gated l2002l l2003l IMivazaki et al.ll2006[ ). Such differences are however 
small, and basically irrelevant compared to the uncertainty introduced by the unexplained strain scaling parame- 
ter 7c that is generally needed to bring the ISHSM into register with experimental data. We refer the reader to 
Rcf. ( Fuchs and Gates 2003 ) for a detailed discu ssion of flow curves from the original ISHSM and the F^^'-modcl, and 



to Refs. (Brader et al 



2005 . [Zausch et aLll2008f ) for additional explicit results. 



VI. CONCLUSIONS AND OUTLOOK 



In our ITT approach to sheared dense suspensions, properties of the stationary state are approximated by following 
the transient structural rearrangements encoded in the transient density correlator <I'q(i) = {Qq^^ * ^?q(^))/-^*S'g = 
{Sq £'q(i) (0) /NSq. Gloscd uou-Markoviau equations of motion for ^q(t), obtained after mode coupling approximations 
(and determined by the equilibrium structure factor, the density and a short time diffusion coefficient) need to be 
solved, while stationary averages under shear are obtained from time- and wavevector integrals containing $q(i). 
Shear advection of density fluctuations accelerates the loss of memory. However, this effect competes with increased 
nonlinear coupling of wavevector modes which models the cage-effect in dense fluids, and which leads to increased 
memory effects and relaxation times. 

ITT proposes a scenario for the transition between a shear-thinning viscoclastic fluid and a yielding/shear-melted 
glass, which captures many features observed in dense colloidal dispersions. The physics of the glass transition has 
thus been brought to bear on addressing the nonlinear rheology of dense dispersions. ITT generalizes the concept 
of an ideal glass transition in (quiescent) MGT to driven (sheared ) systems. It gives stationary averages, correlation 
and response functions (susceptibilities) ([Kriiger and Fuchsl |2008[ ). The essence of the transition between the two 
nonequilibrium states (the elastically distorted glass on the one hand, and the shear-melted state on the other) can 
be captured in schematic models. 

Many open questions still remain in the development and application of ITT-MGT. In the formulation of the 
theory, the neglect of the stress- induced couplings I](i) and E(t) in Eq. ([7^ remains unsatisfactory. We expect it 
is one cause for the theory underestimating the effects of shear, necessitating the introduction of the strain rescaling 
parameter 7c to bring theory and experiment into register. Another contributor to this mismatch is presumably the 
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isotropic approximation underlying the ISHSM; and in Figl2] a third comes from the fitting of incoherent stationary 
correlators to coherent transient ones. Together these factors require a disturbingly small value of 7c to be used, and 
this aspect of the MCT/ITT approach is not yet fully under control. The neglect of the stress couplings leads to 
another aspect of the theory which should be scrutinized further. All yield stress values (the finite stresses remaining 
in the limit 7 ^ in the yielding glassy states) are determined by the equilibrium structure alone. Hence they 
cannot be affected by e.g. hydrodynamic modifications to the bare diffusion constant. The stress-induced couplings 
T,{t) and could change this, as seems required t o explain recent simulations of fluids comprising overdamped 



Newtonian particles under shear (jHorbach and Zausch 20 08). Differences bet ween the stationary and the transient 



density correlators observed in the same simulation study ( Zausch et al.ll2008l) . also are not contained in the present 
description. Recent ly, however, an improved approximation within ITT has been suggested, which comes closer to 



tnproved 

the measured data (jKriiger and Fuchj |2008[ ) . and should in future be used to reconsider the stationary correlators 



from confocal microscopy and simulation discussed in Section [VB I 

Let us finish in pointing to the recent extension of the ITT-MCT approach to arbitrary incompressible and ho- 
mogeneous flows wh ich has led to a general constitutive equation for colloidal dispersions close to a glass transition 
( Brader et al.l[2008h . Consequences of this theory are being worked out. 
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Note added: During completion of this manuscript, we became aware of the work by[Ch ong and KimI (Hooi), who 



worked out an ITT-MCT approach for simple fluids under shear as studied in simulations using the SLLOD equations. 
A detailed comparison of the two approaches may prove very rewarding. 

APPENDIX A: APPENDIX A 

Appendix A contains some technical material used in the text. 

The SO ilj(t) introduced in Eq. (j52p is not perpendicular to density fluctuations and can be analysed using 

^snU ^ gSott (^1 + 7 ^* dt' e-^*' a^y e*"'*'^ = e^* (1 + S](0) . (Al) 

Here, the difference Sfl'^ — = <Jxy between the two advection operators becomes important, which is the shear 
element of the microscopic stress tensor. Its time integral is denoted as Y,{t) in Eq. (jM]) . The SO f^J(i), thus can be 
decomposed as given in Eq. (j53p . 

The time ordered exponential u{t,t') introduced in Eq. ([551) . which satisfies dtu{t,t') = u{t,t')a{t) with u{t,t) ~ 1, 
and where a{t) is an operator that docs not commute with itself for different times, is given by (for t > t'): 

u{t,t') = e^^''^'"^'^ = 1 + ^ [ dsi [ ds2 ... [ ds„ a(s„)a(s„_i)...a(si) . (A2) 

,1=1"'*' 

It obeys: u{t,t') = u{t" ,t') u{t,t") for any t > t" > t', dt'u{t,t') = -a{t')u{t,t'), and u{t,t')-^ = v{t,t') = e~^'''^'''^''\ 
where e+ is built w ith the reverse order o f operators compare d to E q. (|A2[) . and solves dtv{t,t') ~ —a{t) v{t,t') with 
v{t,t) = 1; see e.g. (jKawasaki and Gunton 1973 , van K ampcn. ll98l[ ). 

The irreducible dynamics appearing in the final formally exact equation of motion Eq. ((551) ; can be rewritten using 
Eq. (|7T1) . Standard manipulations starting from dtUi{t,t') = Ui{t,t') fj|(<) lead to 

u,it,t') = u^it,t') + f dt' u,it',o) nl{t')Qit') u^it,t') 







where the newly defined time evolution operator 

U^{t,t') = cxp_ { ds n^is) Q{s)} ^ exp_ { ds e^"'^ Q{s) nl er^^'' Q{s)} 
satisfies P U^{t, t') ^ P = uf^{t, t') P. This leads to Eqs. (|73l74p under the approximation of Eq. 
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APPENDIX B: APPENDIX B 



The ITT approach introduced bv lCates et alj ( 2004[ ). lFuchs and CatesI ( 20021 [20051) . consisted of equations of motion 
shghtly different from the ones derived here. Nevertheless, it led to the identical nonlinear stability equation Eq. ((M)) 
and thus to the identical universal scenario of a shea r thinning fl uid and a yielding glass, wh ich is also the basis and 
content of the schematic models bv lFuchs and CatesI (|2002L l2003l) . The ISHSM presented bv lFuchs and CatesI (|2002l 
I2003D differed from the one summarized in Appendix C in quantitative details, but by an amount swamped by the 
unknown parameter 7c entering because of isotropic averaging. Thus detailed quantitative comparisons between model 
variants are not warranted. Nonetheless, in this Appendix, the differences between the microscopic ITT formulations 
of the different papers will be listed to establish an unambiguous formulation. 

A first difference concerns the definition of the transient correlator, which in version (I) of ITT (jCates et al.ll2004l 
iFuchs and c"ate^l2002L l20?)5h carried the wavevector advection on the left hand side of the average 



NSa 



(£'q(-t) e filq) 



Sn 



where Eq. ([28)) was used. This definition, and the present one in Eq. ([28]) obey translational invariance and at first 
sight seem equivalent. Obviously, the sign of 7 is arbitrary, so that the difference between forward and reverse 
advection appears irrelevant. Yet, the version (I) correlator obeys for short times 



<i>m(i^o) 



1 



Sq ^ q 



Sn 



Sn 



qxqy 



which shows that the initial decay rate in ITT- (I) is affected by the stretching of equilibrium density fluctuations 
contained in Sg(^t) . This subtle difference is the reason for using the present dcflnition of the advected wavevector 
in Eq. (|15p . This complication could be got rid off in ITT-(I) by performing the limit of 7 — *■ 0, accompanied by 
t — > 00 keeping the strain ^7 = const.. However, the present formulation, where 7 always appears in the fashion of 
an accumulated strain -jt, appears simpler. Moreover, it has not been worked out yet, how the second difference in 
the ITT equations of motion, to be discussed next, can be eliminated from version (I). Therefore, we suggest to use 
deflnition Eq. (|28p for the transient correlator from now on. 

An important difference to ITT- (I) concerns our present handling of the linear coupling of density fluctuations to 

themselves. It leads to an 'initial decay rate' Tq{t) that becomes time-dependent under shear. In ITT-(I), Tq\t) 
could exhibit different signs depending on the direction of q. In numerical solutions of ITT-(I), the appearance of 

1/Tq\t) as prefactor in the memory kernel destabilized the numerical algorithm at long times, restricting the density 
and shear rate window where ITT- (I) could be applied (.Hc nrich 2007b ) . The stability of the equations of motion 
derived bv lCates et al.l ( 20041 ) . iFuchs and CatesI ( 20021 . |2005| ) could not be assured. While stable equations might be 
possible with the old definition of the transient correlator, the new definition in Eq. (pS)) easily led to the convenient 
time evolution operator in Eq. (j48p . Again, it appears not to be the only choice giving stable relaxations, but serves 
well at least as long as the accumulated elastic energies E(t) are neglected. Introduction of the advection operator 
jri''', working symmetrically with collective densities, and the special choice of projection in Eq. (|5T|) overcomes the 
problem of negative rq(t), and leads to the vertices in the memory function of Eq. (|ST|) . which assure a positive 
instantaneous friction. Therefore, we suggest to use the approxi mation Eg. (1811) for the memory kernel from now on. 
Quantitatively, but not qualitatively, the results on step-strain bv lBrader et al. need to be recalculated, because 



there the details of the ITT- ( I) ver t ex ma ttered. All other results presented in Refs. ( Brader et al.l 2007 



2004', 'Fuchs and Catcs"2002', '2003', '2005') remain unaffected, as do all results in Refs. (jCrassous et al.l 



Cates et aLl 
20061 . 12OO8I 



Fuchs and Cates 2002 , ,2003. .Hainal and Fuchs 2008 . Varnik and Henrich|[2006[) based on schematic models. However, 
as already mentioned, the result given for the isotropic term in the distorted structure factor bv lHenrich et al. I (I2OO7D 
was inaccurate as it assumes consistency at large q which ITT-MCT in fact docs not achieve. 



APPENDIX C: APPENDIX C 



For completeness, the numerically solved ISHSM is summarized. It consists of Eq. ([5^ with the assumption of 
isotropic correlators, ^q{t). The initial decay rate is approximated without shear; = q^Dg/Sq, where Dg is the 
density-dependent short time diffusion coefficient. The memory function also is taken as isotropic and modeled close 
to the unsheared situation 

"^^W ^ ^E^qJW Mt)<^\^-i.\it) , (CI) 

k 
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FIG. 5: Integrands for the parameter c^'^^ in the stabihty equation (|84|) : the one labeled c^°'=° resuhs for the memory function 
in Eq. (|8TJ, the one labeled c'Shsm ISHSM approximation Eq. (fCT|) setting 7^ = 0.84 to recover c'"^' = 0.64 in both 

cases. For comparison the integrand leading to the exponent parameter A is included as dotted line. 



with 



1 SqSk.Sp 



q • k Cfe(4) + q • p Cp(t) [q • k Cfc + q • p Cp] 



(C2) 



where p = q— k, and the length of the advected wavevectors is approximated by k{t) = k{l + {t'^ HcYY^"^ and p{t) = 
p(l + (^7/70)^)^^^ • Note, that the memory function thus only depends on one time. We use the phenomenological 
adjustment factor 7^ in order to correct for the underestimate of the effect of shearing in the ISHSM. 
The expression for the potential part of the transverse stress may be simplified to 



dt dk 



607r2 



while the (anisotropic) structure factor is approximated as 



(C3) 



^ ^aniso 



(7) 
(7) 



7c JO 



dtj^^mit) for 7r«l 



dt '^'^ 



m 



for 1 



(C4) 



In the last two equations, Eqs. (|C3IC4p . the parameter 7c is taken to be 7c = -s/S, as would follow from isotropic 
averaging of k(t). Note that in the glass, the distorted structure factor becomes isotropic in ISHSM, because for 
7T 3> 1, the isotropization approximation should be performed on Eq. (|88p . which leads to the second line in Eq. 
(IC4l) 



For the numerical solution of the ISHSM the algorithm bv lFuchs et al] (Il99l|l was used for ha rd spheres with radius 
R. Their structure factor Sk is taken from the Percus-Yevick approximation ( Russel et al. 19891 ) and d epends only on 
the packing fraction (f>. The wavevector integrals were discretized according to iFranosch et al. [ (|l997f ) with M = 100 
wavevectors from fcmin = 0.1/i? up to /cmax = 19.9/i? with separation Afc = 0.2/R. Time was discretized with 
initial step-width dt = 2 x 10~''i?^/£'o, which was doubled each time after 400 steps. We find that the model's glass 
transition lies at (pc = 0.51591. Thus e = {4> — 4>c)/4>c (where e ~ 1.54e), and 7 are the only two control parameters 



determining the rhcology. The exponent parameter becomes A 
still change somewhat if the discretization is made finer. 



0.735 and c^^^ w 0.45/7c ; note that these values 
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The ISHSM value for c^^^ gives a rough estimate of the error in ISHSM caused by using isotropically averaged 
vertices. Without this approximation, c^^-* ~ 0.64 holds for the present discretization in an anisotropic calculation. 
This leads to the conclusion that 7c ~ 0.84 should approximately correct for the error when isotropically averaging. 
Yet, as discussed in Section |VB] appreciably smaller values 7c < 0.1 are required to fit the measured data. More 
information on the spatial distribution of the erro r in the ISHSM can be gained from considering the integrands in 
the expressions leading to c^^') = J dqi^^ derived bv lFuchs and Gated (|2003). They are shown in Fig. [5] for the ISHSM 
using the value 7c = 0.84, and for the anisotropic calculation. The ISHSM not only errs in magnitude, (e.g. in the 
prefactor of the shear induced relaxation time) but also qualitatively misses the negative contributions. The curve 
with anisotropic shear-advection taken properly into account indicates that the particles rearrange locally during the 
initial part of the yield process, as contributions from the range 2 < qR < 20 dominate. 
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